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The  present  dissertation  is  motivated  by  a  desire  to 
have  a  more  precise  knowledge  of  asymptotic  approximation 
error  than  that  given  by  best  order  of  approximation.  It 
owes  its  inspiration  to  a  paper  by  G.  Birkhoff  and  A.  Priver 
concerning  error  bounds  for  derivatives  of  Hermite 
interpolation  and  a  paper  of  C.  A.  Hall  and  W.  W.  Meyer 
concerning   error    bounds    for    cubic    splines. 

In  Chapter  One  we  consider  well  known  results 
concerning  interpolation,  polynomial  approximation  and 
error  analysis  of  spline  approximation.  The  results  given 
here  are  meant  to  provide  a  context  for  the  theorems  given 
in  later  chapters.  In  Chapters  Two  and  Three  we  consider 
the  problem  of  best  error  bounds  for  derivatives  in  two 
point   Birkhoff    interpolation   problems. 


Chapter  Four  presents  the  problems  of  existence, 
uniqueness,  explicit  representation,  and  the  problem  of 
convergence  for  fourth  degree  splines.  Moreover  we  also 
consider  the  problem  of  optimal  pointwise  error  bounds  for 
functions  f  e  C^^^  [0,1].  In  Chapter  Five  our  main  object 
is  to  sharpen  the  error  bounds  obtained  earlier  by  Marsden 
concerning  quadratic  spline  interpolation.  By  doing  so  we 
obtain  in  some  special  cases  error  bounds  that  are  in  fact 
optimal. 


CHAPTER  ONE 
INTRODUCTION 


The  purpose  of  this  chapter  is  to  provide  a  context 
for  the  results  derived  in  succeeding  chapters.  In  order 
to  show  some  of  the  important  achievements  in 
approximation  by  polynomials,  we  discuss  briefly  the 
Lagrange  and  Hermite-Fe jfer  interpolations,  which  match  a 
given  function  at  any  finite  number  of  distinct  points. 
After  exploring  the  question  of  computational  stability  of 
a  given  interpolation,  we  discuss  in  some  detail  the 
problem  of  best  order  of  approximation  by  polynomials  as 
initiated  by  S.  N.  Bernstein  [1912],  D.  Jackson  [1930], 
and  A.  Zygmund  [1968]. 

In  contrast  to  high  order  approximation  by  a  single 
polynomial,  we  next  consider  in  great  detail  the  problem 
of  approximating  a  given  function  f(x)  defined  on  [a,b] 
by  the  interpolatory  piecewise  polynomials  known  as 
splines.  Special  attention  is  given  to  the  problem  of 
approximating  by  piecewise  cubic  and  piecewise  parabolic 
splines.  The  study  of  these  splines  motivates  us  to  also 
study  two  point  Hermite  and  Birkhoff  interpolations. 


Lagrange   and   Hermite-Fe j^r    Interpolation 
Let  X  denote  an  infinite  triangular  matrix  with  all 
entries    in    [-1,    1] 

^0,0 

(1.1.1)  X:         Xo;i      xi^i 

^0,n   ^l,n   •     •         ^n,n- 
We  denote  by  Lj^[f,x;X]  the  Lagrange  polynomial  of 

interpolation  of  degree  _<  n   which  coincides  with  f(x)  in 

the  nodes  Xj^^^  (k  =  0 ,  1 ,  .  .  ,  n)  .   Then 

n 

(1.1.2)  Lj^  [f,x;X]  =  I      f{X]^^)l]^^{x) 

k  =  0 

where 

(1.1.3)  lv„(x)  =     '^n(^)     , 


•kn' 


(X  -  Xvn)   t^n'  (Xvn) 


kn'  "^n  ^^kn' 


n 
coj^(x)  =   n   (X  -  x^^)     . 
k=0 

It  is  known  from  the  results  of  G.  Faber  and  S.  N. 
Bernstein  that  no  matrix  X  is  effective  for  the  whole 
class  C  of  functions  continuous  in  [-1,  1].  Bernstein 
showed  that  for  every  X,  there  exists  a  function  fo(x)  and 
a  point  Xq  in  C[-l,l]  such  that 
(1.1.4)    lim^_>  ^  Lj^[fo,Xo;X]  =  CO   . 

L.  Fejfer  [1916]  showed  that  if  instead  of  Lagrange 
interpolation,  we  consider  the  Hermite-Fe jfer  interpolation 
polynomials,  the  situation  changes.  The  Hermite-Fe jfer 
polynomials  H^^^ -|_  [  f  ,x  ,X]  are  of  degree  <_2n  +  1  and  are 
uniquely  determined  by 


(1.1.5)  H^^i[f,x^^;X]     =f(xj,^),       H^41)[f,Xj,^;X]     =       &^^ 
where    6]^^    are    arbitrary    real    numbers,       k    =    0,    1,    .    .    n. 
The    explicit   form   of   Hj^^-|_  [f  ,x;X]    is   given  by 

(1.1.6)  H^^i[f,x;X]     =      I       f(xj^^)    hj,^(x)     +      I       j,^        j,^(x) 

k=0  k=0 

where 

(1.1.7)  hj^n^^^    =    [1    -  '^n'  '  <^kn^  ^^   "   ^kn^  }         l^n^  ^^^ 

ui^Mxj^^) 

=  •   ^kn<^)    Ikn^(^) 
and 

(1.1.8)  n^^(x)     =     (X    -    X^^)     l^^^(K)      . 

Fejfer  brought  out  the  importance  of  Hermite  interpo- 
lation by  introducing  the  concept  of  "strongly  normal" 
point  systems.  To  each  set  of  n  +  1  distinct  points  Xq, 
Xi,  .  .  ,  X  ,  Fej^r  associates  a  set  of  n  +  1  points  Xq, 
X-,,  .  .  ,X^  which  are  the  zeros  of  the  linear  functions 
j^(x).  The  points  Xq,  X-|^,  .  .  .  ,  X^^  are  said  to  be  the 
conjugate  point  system  of  Xq,  x-j_,  .  .  ,  x^.  A  system  of 
points  Xq,  x-^,  •  .  ,  Xj^  is  called  strongly  normal  if  the 
conjugate  point  system  lies  inside  [-1,  1].  For  example, 
the  zeros  of  the  Tchebycheff  polynomial  T^(x)  =  cosnG, 
cosG  =  X  form  a  strongly  normal  point  system.  Fejir 
proved  (using  these  ideas)  that  Hermite-Fe jfer  interpola- 
tion polynomials  based  on  strongly  normal  point  systems 
(and  under  certain  conditions  on  5]^^^)  converge  uniformly 
to  f (x)  on  [-1,  1] . 


Optimal  Error  Bounds  for  Two  Point  Hermite  Interpolation 
In  order  to  motivate  the  present  day  work  on  error 
bounds,  we  first  consider  the  classic  error  bound  of 
Cauchy.  Let  us  consider  once  more  the  interpolation 
formula  of  Lagrange.  Let  f(x)  e  C[a,b]  and  consider  the 
Lagrange  interpolation  polynomial 

^n  ff'^J  =   I  f(^kn)  Ikn(^)  • 
k=0 

Next  we  set 

(1.2.1)  e(x)  =  f(x)  -  L^  [f,x]  . 

In  the  case  f  (x)  is  itself  a  polynomial  of  degree 
<_  n,  then  it  is  easy  to  see  from  the  uniqueness  of  the 
Lagrange  interpolation  polynomial  that  e(x)  =  0.  Thus  it 
is  of  interest  to  study  what  can  be  said  about  e(x)  if 
f(x)  is  a  given  smooth  function  other  than  a  polynomial  of 
degree  <_  n.  The  following  theorem  gives  the  most  widely 
known  error  bound. 

Theorem  1.1  (Cauchy).   Let  f(x)  e  C[a,b]  and  suppose 
that  f^'^'(x)  exists  at  each  point  of  [a,b].   Let  Lj^[f,x] 
be  the  element  of  the  class  of  polynomials  of  degree 
j<  n  -  1  that  satisfies  the  equation 

(1.2.2)  L^[f,x^^]  =  f(Xi^)   i  =  0,  1,  .  .  ,  n. 
Then  for  any  x  in  [a,b],  the  error 

e(x)  =  f(x)  -  Lj^[f,x] 
has  the  value 

(1.2.3)  e(x)  =    i^^(x}    f  ('^^^M  ?  )/(n  +  l)  !  , 

where  ?  is  a  point  of  [a,b]  that  depends  on  x  and 


n 
i  =  0 
An   immediate   consequence   of    (1.2.3)    is    the    inequality 
(1.2.4)  |e(x)  I    <    |w^(x)  I     I  |f  (^-^l)  I  |/(n+l)  ! 

where  | |  | |  denotes  the  supremum  norm  on  [a,b].  If  we  set 
f(x)  =  oj  (x)  ,  we  see  that  (1.2.4)  becomes  an  equality. 
Thus  the  right  hand  side  cannot  be  made  smaller.  We 
therefore    say   that    (1.2.4)    is    an   optimal    bound. 

The  Equations  (1.2.3)  and  (1.2.4)  have  been 
extensively  studied.  For  instance,  the  study  of 
minimizing  ll'^nN  ^^'^  ^°  Tchebychev's  system  of 
orthogonal  polynomials.  For  a  good  discussion  of  some  of 
the  elementary  analysis  associated  with  this  error  bound, 
see    Powell     [1981]. 

In  contrast  to  the  precise  and  beautiful  pointwise 
Cauchy  bound,  very  little  has  been  known  about  precise 
polynomial  derivative  errors.  Denoting  e(x)  as  the  Cauchy 
remainder  for  Lagrange  polynomial  interpolation,  we 
consider  the  role  played  by  the  term  f  *^"*'-'-^  (^) .  If  f  e  Pj^ 
(the  class  of  polynomials  of  degree  <_  n)  ,  the  remainder 
vanishes  identically.  For  a  fixed  x,  we  may  consider  the 
remainder 

e^(x)    =    f(x)    -   Lj^[f,x] 
as  a  process  which  annihilates  all   elements   of  P    .      We 
may   now   formulate   the   following   theorem   of   Peano    [1913]. 


Z^££i^IE  ^'  ^     (Peano).   Let  L  be  a  continuous  linear 

functional   such  that  L(p)  =  0  for  all  p  6  Pj^.   Then  for 

all  f  6  c^"^^)  [a,b]  , 

b 
(1.2.4)   L(f)  =   /  f^^  +  ^Mt)  K(t)  dt 
a 

where 

K(t)  =  [  Lj^[(x  -  t)^]     3  /  n! 

and 


(X  -  t)^ 


=    (X  -  t)"     for  X  >_  t 
=  0  for  X  <  t  . 


The  notation  L  [(x-t)"]  means  that  the  functional  L 
is  applied  to  (x  -  t) ^  considered  as  a  function  of  x. 
For  a  detailed  study  of  the  Peano  theorem  we  refer  to  P. 
J.  Davis  [1975]  and  to  A.  Sard  [1963].   We  next  turn  to  an 
application   of  the  Peano  theorem  to  derive  pointwise 
optimal  derivative  error  bounds. 

Let  u(x)  e  C^'^^[0,  h]  be  given;   let  V3  (x)  be  the 
unique  Hermite  interpolation  polynomial  of  degree  <_  3 
satisfying 
(1.2.5)    V3(0)=u(0)  ,  V3(h)  , 

v'3(0)=u'(0)  ,        v'3(h)=u'(h)  . 

Ciarlet,  Schultz  and  Varga  [1967]  obtained  a 
pointwise  error  bound  for  e  (x)  =  V3(x)  -  u(x)  and  its 
derivatives    in   terms   of 

U    =    maxo<x<hh^'^^  (^)  I     • 


Their  bounds  are 

(1.2.6)  |e^^Mx)|  <    h^  [x(h-x)]^~^  U   k  =  0,  1,  2. 

k!  (4  -  2k)  ! 

For  k  =  0,  (1.2.6)  is  best  possible,  since  equality 
holds  for  u(x)  =  x'' (h-x)  ,  whose  Hermite  interpolation 
polynomial    is    v=0. 

G.     Birkhoff     and    A.     Priver     [1967]     obtained    the 

I     (kl  I 

following   optimal    error   bounds    on    the    derivative    |e^    '  (x)  | 

in  terms  of  U. 

Theorem  1.3  (Birkhoff  and  Priver).  Let  u(x)  G  c'^CO,!]. 

Then  we  have  (h  =  1) 

(1.2.7)  |e' (x) |/U  <  [  x(x-l) (2x-l)  ]  /  12 

for  0  <  X  <  1/3  , 
<_  [  16x^  -  lOSx^  +  197x  -  162 

+  66/x  -  13/x2  +  l/x^  ]  /  96 

for  1/3  <  X  <  1/2  . 

(1.2.8)  |e' ' (x) |/U  <     [    48x^  +  42x^  -  lOOx^ 

+  54x2  -  12x  +  1  ]  /  2(l-x)^ 
for  0  _<  X  j<  1/3  , 
<  [  -6(x-l/2)2  +  1/2  ]  /  12 

for  1/3  <  X  <^  2/3  . 

(1.2.9)  le' * ' (x) |/U  <  -(x-1/2)^  +  3(x-l/2)2/2  +  3/16 

for  0  _<  X  £  1  . 
For  1/2  <_  X  <_  1  the  bounds  of  e^^^  (x)  are  given  by 

(1.2.10)  e^^Mx)  =  e^^^  (1-x)   ,     k  =  0  ,  1,  2,  3. 
Further,  from  Birkhoff  &  Priver,  the  uniform  error 

bounds  are  given  by 


e^^^  (X)  I  <  aj.  U  ,   r  =  1,  2,  3, 


4^4! 
(1.2.11)  a^L  =   (/3)/216 

32  =  1/12 

33  =  1/2  . 

The  proof  of  the  above  theorem  is  based  on  the  Peano 
kernel  theorem.  It  gives  a  general  and  highly  useful 
method  for  expressing  the  errors  of  approximations  in 
terms  of  derivatives  of  the  underlying  functions  of  the 
approximation.  For  a  computer  routine  which  gives 
polynomial  error  bounds  by  numerical  quadrature  of  the 
Peano  kernel,  see  Howell  and  Diaa  [1986].  Stroud  [1974] 
gives  a  readable  account  of  some  other  applications. 

Birkhoff  Interpolation 

We  have  just  observed  that  in  problems  of  Hermite 
interpolation,  function  values  and  consecutive  derivatives 
are  prescribed  for  given  points.  In  1906,  G.  D.  Birkhoff 
considered  those  interpolation  problems  in  which  the 
consecutive  derivative  requirement  can  be  dropped.  This 
more  general  kind  of  interpolation  is  now  referred  to  as 
the  Birkhoff  (or  the  lacunary)  interpolation  problem (s). 

The  Birkhoff  interpolation  problem  differs  from  the 
more  familiar  Lagrange  and  Hermite  interpolation  in  both 
its  problems  and  its  methods.  For  example,  Lagrange  and 
Hermite  interpolation  problems  are  always  uniquely 


solvable  for  every  choice  of  nodes,  but  a  given  Birkhoff 
interpolation  may  not  give  a  unique  solution. 

More  formally,  given  n  +  1  integer  pairs  (i,k) 
corresponding  to  n  +  1  real  numbers  c^  y^,  and  m  distinct 
real  numbers  x^^,  i  =  1 ,  2 ,  ,  ,  m  <_  n  +  1 ,  a  given  problem 
of  polynomial  interpolation  is  to  satisy  the  n  +  1 
equations 

(1.3.1)    Pn^^Nx^)  =  y^^j, 

with  a  polynomial  P^  of  degree  at  most  n.   (We  are  using 
the  convention  that  P^^^^  (x)  =  Pj^(x).) 

If  for  each  i,  the  orders  k  of  the  derivatives  in 
(1.3.1)  form  an  unbroken  sequence  k  =  0,  1,  .  .  ,k^,  then 
the  interpolation  polynomial  always  exists,  is  unique,  and 
can  be  given  by  an  explicit  formula.  If  some  of  the 
sequences  are  broken,  we  have  Birkhoff  interpolation.  As 
remarked  by  Professor  Lorentz  [1983],  the  two  cases  are  as 
different  as,  let  us  say,  the  theory  of  linear  and 
nonlinear  differential  equations. 

Pairs  (i,k)  which  appear  in  (1.3.1)  are  most  easily 
described  by  means  of  the  interpolation  or  incidence 
matrix  E.  If  Pj^^'^'(Xj_)  is  specified  in  (1.3.1),  we  put  a 
"1"  in  the  i  +  lst  column  and  kth  row  of  E.  If  P^^*^'  (xj_)  is 
not  specified  in  (1.3.1),  then  a  "0"  appears  in  the  i  +  lst 
column  and  kth  row.  Each  of  the  m  rows  of  E  has  a  non- 
zero entry.  An  incidence  matrix  E  and  a  pointset  X,  which 
lists  the  points  Xj_,  specify  a  Birkhoff  interpolation 
problem  of  the  type  of  (1.3.1).   For  a  given  E  and  X,  the 


10 


unique  existence  of  an  interpolation  polynomial  of  degree 
n  +  1  is  equivalent  to  the  invertibility  of  the  system  of 
equations  given  by  (1.3.1),  or  equivalently  to  the  inver- 
tibility of  a  matrix  V  which  we  will  refer  to  as  a 
generalized  Vandermonde  matrix  V.  For  Lagrange 
interpolation  of  the  points  Xj_,  i  =  1 ,  2 ,  .  .  ,  n  +  1 , 
the   Vandermonde   V   is   given   as 

1  1         .       .       .       .  1 


1.3.2)         V 


x-L         X2       .       .  Xj^+2. 


„   n  n  V      i^ 


'1         ^2       •       •  ^n+1    • 

Inversion  of   the   Vandermonde   gives    the   coefficients   of   the 

fundamental    functions    lj,j^(x)    of  Lagrange   interpolation. 
As   Lagrange   interpolations   are   always   unique,    it    follows 
that   Vandermonde    matrices    are    invertible. 

For  a  given  system  (1.3.1),  it  is  not  hard  to 
construct  an  analagous  matrix  to  (1.3.2),  which  we  will 
refer  to  as  the  generalized  Vandermonde.  Just  as 
inverting  the  Vandermonde  matrix  gives  the  fundamental 
functions  of  Lagrange  interpolation,  inverting  the  genera- 
lized Vandermonde  gives  a  convenient  form  for  representing 
a  Birkhoff  interpolation.  The  Vandermonde  and  its 
counterpart  for  Birkhoff  interpolation  are  examples  of 
Gram  matrices,  of  which  a  good  account  is  to  be  found  in 
Davis     [1975]. 

Though   invertible,    the  Vandermonde   matrices   are   known 
to  be  extremely   ill-conditioned   for   real-valued 
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interpolation.  Many  of  the  generalized  Vandermonde 
matrices  associated  with  Birkhoff  interpolation  processes 
are  much  better  conditioned,  illustrating  an  advantage  of 
Birkhoff  interpolation  over  the  more  traditional  Lagrange 
interpolation.  To  make  this  point  more  explicit,  we 
define  "condition"  of  a  matrix. 

For  a  given  norm  |  |  |  |,  and  invertible  matrix  M,  we 
define  the  condition  cond(M)  of  the  matrix  M  by 
(1.3.3)    cond(M)  =  | |M| J  | |m"^| |  . 

If  we  rescale  the  Birkhoff  interpolation  problem 
specified  by  E  and  X  to  the  unit  interval,  we  can  define 
the  condition  of  an  interpolation  as  the  condition  of  the 
associated  generalized  Vandermonde.  In  the  L2  norm  for 
eleven  equally  spaced  points,  the  condition  number  of 
Lagrangian  interpolation  is  on  the  order  of  a  million.  On 
the  other  hand,  Lagrangian  interpolation  on  eleven  equally 
spaced  complex  roots  of  unity  has  L2  condition  number  one, 
as  does  the  eleven  term  MacLaurin  expansion. 

Computationally  speaking,  the  inverse  of  the 
condition  number  of  a  matrix  M  is  the  norm  distance  of  M 
from  a  singular  matrix  (See  Golub  and  Van  Loan  [1983]). 
For  example,  the  Vandermonde  for  Lagrange  interpolation  of 
eleven  points  on  the  unit  interval  is  thus  seen  to  be  a 
norm  distance  of  only  one-millionth  from  being  singular. 
Not  only  is  the  il 1-conditionedness  of  the  Vandermonde 
troublesome  in  determining  the  coefficients  of  the 
fundamental  functions,  but  it  also  causes  problems  of 
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round-off  error  in  evaluating  a  polynomial  by  use  of  the 
fundamental  functions.  For  these  reasons,  it  is  very  much 
preferable  to  use  a  well-conditioned  interpolation. 

The  MacLaurin  expansion,  having  diagonal  generalized 
Vandermonde,  is  as  well-conditioned  as  is  possible. 
Another  particularly  well-conditioned  interpolation  is  the 
Lidstone  interpolation. 

A  Lidstone  polynomial  is  a  truncation  of  a  Lidstone 
series.  In  turn,  a  Lidstone  series  is  a  generalization  of 
a  Taylor  series  which  approximates  a  given  function  in  the 
neighborhood  of  two  points  instead  of  one.  Such  series 
have  been  studied  by  G.  J.  Lidstone  [1930] ,  by  Widder 
[1942],  by  Whittaker  [1934]  and  by  others.  More 
precisely,  the  series  has  the  form 

(1.3.3)  f(x)  =  f(l)Ao(x)  +  f(0)AQ{l-x)  +  f'(l)Ai(x)  + 

f ' '  (O)A^(l-x)    +    .    .    . 
where   Aj^(x)    is   a   polynomial   of   degree    2n    +    1    defined   by 
the   relations 

Aj^(x)    =   X 

(1.3.4)  A^' ' (X)     =    An-l(^) 

A^(0)  =  A^(l)  =0,  n  =  1,  2,  .  .  . 
Thus  it  is  clear  that  the  sum  of  an  even  number  of 
terms  of  the  series  (1.3.3)  is  a  polynomial  which  coin- 
cides with  f(x)  at  X  =  0  and  at  x  =  1.  Moreover,  each 
even  derivative  of  the  polynomial  coincides  with  the 
corresponding  derivative  of  f(x)  at  those  points. 
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Polynomial  Approximation 

Weierstrass  first  enunciated  the  theorem  that  an 
arbitrary  continuous  function  can  be  approximately 
represented  by  a  polynomial  with  any  degree  of  accuracy. 
We  may  express  this  theorem  in  the  following  form. 

If  f(x)  is  a  given  function,  continuous  for 
a  <_  X  <_  b,  and  if   e   is  a  given  positive  quantity,  it 
is  always  possible  to  define  a  polynomial  P(x)  such  that 

(1.4.1)  |f  (X)  -  P(x)  I  <  e 
for  all  a  <_  X  <_  b. 

It  is  readily  seen  that  the  number  of  terms  required 
to  yield  a  specified  degree  of  approximation,  or  under  the 
converse  aspect,  the  degree  of  approximation  attainable 
with  a  specified  number  of  terms,  is  related  to  the 
properties  of  continuity  of  f(x).  Naturally  this  has  led 
to  many  interesting  developments  in  the  theory  of  degree 
of  approximation  of  continuous  functions  by  polynomials  to 
which   we    turn   to   describe. 

A  first  important  step  in  building  this  theory  was 
made  by  D.  Jackson  [1930].  Let  f  G  C[-l,l].  Suppose  that 
we  define  the  best  approximation  of  f  by  polynomials  of 
degree   n   by 

(1.4.2)  Ej^(f)     =    inf    I  |f    -    Pj^l  I 

where  P  ranges  over  all  algebraic  polynomials  of  degree  n 
and  1  I  f  I  I  =  max  |  f  (x)  |  ,  a  <_  x  <_  b.  Jackson  considered  the 
problem  of  estimating  E^(f).  To  describe  his  results  we 
need  the  following  definition. 
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Definition  1.1   If  f  e   C[a,b],  then  the  modulus  of 
continuity  of  f  is  a  function   (f,h)  such  that 

(1.4.3)  (f,h)  =  sup|^_y|<h;  x,y  6  [a,b]|f(^)  "  ^^^H  ' 
Now  Jackson's  theorems  may  be  easily  stated. 

Theorem  1.4  (Jackson).  Let  f  be  continuous  on  [-1,1]. 
There  is  a  positive  constant  A  such  that 

(1.4.4)  Ej^(f)  <  A  co(f,l/n)  ,  n  =  1,  2,  .  . 
where  A  is  independent  of  f. 

An  important  corollary  of  Theorem  1.3  deserves  to  be 
mentioned.  Let  Lipa  [■_]_  ^^  j  (M)  (or  simply  Lipa )  be  the 
class  of  functions  f  in  C[-l,l]  such  that 

If  (X)  -  f  (y)  I  i  M  |x-yi^ 
for  all  X  and  y  in  [-1,1].   It  is  easy  to  see  that 
f  6  Lipctr^  ^■,   (M)   if  and  only  if 

co(f,h)  <_  M  h^'for  all  h  >_  0  . 
We  then  have  the  following  consequence  of  Jackson's 
theorem. 

Corollary  1.5   Let  0  <  a<^  1.   If  f  6  Lipc'[_i^i]  (M),  for 
some  constant  M,  then 

(1.4.5)  Ej^(f)     <   _A_        for   n    =    1,    2,     .     .     . 

n^ 
for   some   positive   constant   A. 

A.       F.       Timan       [1951]      noticed     the      following 

strengthening    of    Jackson's    theorem. 

Theorem   1.6    (Timan).      There  is   a  positive  constant  C  such 

that  if  f  e    C[-l,l]    and  n  is  a  natural  number,   then  there 
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is  a  polynomial  P^^  of  degree  n  such  that 


(1.4.6)  |f(x)    -    Pn(x)|     <      A[    aj(fyi    -   x^)     +   a)(f,l/n2)     ] 

n 
for  all  X  in  the  interval  [-1,1]. 

In  this  result,  in  contrast  to  the  theorem  of 
Jackson,  the  position  of  the  point  x  in  the  interval 
[-1,1]  is  taken  into  consideration  and  it  is  apparent  that 
for  the  polynomial  P^^^^  thus  constructed,  as  |x|  ->  1, 
the  deviation  |f(x)  -  P  (x)  |  is  of  magnitude  a)(f,l/n  ). 

Following  the  important  theorem  of  Timan,  V.  K. 
Dzjadyk  [1956]  proved  the  converse  of  Jackson's  theorem. 
Theorem  1.7  (V.  K.  Dzjadyk).   Let  f  6  C[-l,l].   Suppose 
that  0  <  a  <  1.   Then  there  is  a  constant  B  such  that  to 
each  n  there  corresponds  a  polynomial  Pj^  of  degree  n  such 
that 


(1.4.7)    |f  (x)  -  P_(x)  I  <   B'  [  (  /I  -  x^  )^+    (  1  )^J 

2 
n  n^ 

if  and  only  if  oj  (f  ,h)  £  C  h^  for  some  constant  C. 

From  Jackson's  theorem  we  noticed  that  if  f  G  Lipa, 

then 

E^(f)     <      AM    '      n    =    1,    2    .     .     . 
n^ 
where    A    is    an    absolute    constant.      To    achieve    a    more    rapid 

decrease    to    0    of    Ej^(f),    it    is    necessary    to    assume    more 

smoothness     for     f,     for    example,     that    f    has     several 

continuous    derivatives.       Let    C-'^[-l,l]     ,     r    =    0,     1    ..     . 

denote    the    subset   of   C[-l,l]    consisting   of   those   functions 
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which  possess  r  continuous  derivatives  on  [-1,1],  For 
this  class  of  functions,  Dunham  Jackson  proved  also  the 
following  direct  theorem. 

Theorem  l.B  (D.  Jackson).  If  f  G  c'^^[-l,l],  then 
(1.4.8)    E^(f)  <  Aj,  (1/n)^   a3(f^^^l/n)  ,  n  =  1,  2,  .  .  . 
For  many  important  contributions  we  refer  to  the  work 
of  G.  G.  Lorentz  [1983]  . 

Spline  Approximation 

One  uses  polynomials  for  approximation  because  they 
can  be  evaluated,  differentiated  and  integrated  easily  and 
in  finitely  many  steps  using  just  the  basic  arithmetic 
operations  of  addition,  subtraction  and  multiplication. 
But  there  are  limitations  of  polynomial  approximations. 
For  example,  the  polynomial  interpolant  is  very  sensitive 
to  the  choice  of  interpolation  points.  If  the  function  to 
be  approximated  is  badly  behaved  anywhere  in  the  interval 
of  approximation,  then  the  approximation  is  poor  every- 
where . 

This  global  dependence  on  local  properties  can  be 
avoided  when  using  piecewise  polynomial  approximation. 
Concerning  piecewise  polynomial  approximation,  Professor 
I.  J.  Schoenberg  remarked  that  "polynomials  are  wonderful 
even  after  they  are  cut  into  pieces,  but  the  cutting  must 
be  done  with  care.  One  way  of  doing  the  cutting  leads  to 
the  so-called  spline  functions"  (Schoenberg  [1946], 
p.  46)  . 
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Splines  were  introduced  by  Prof.  Schoenberg  in  1946 
as  a  tool  for  the  approximation  of  functions.  They  tend 
to  be  smoother  than  polynomials  and  to  provide  better 
approximation  of  low  order  derivatives.  Though  v/e  will 
later  use  the  word  spline  in  a  somewhat  broader  context, 
we  first  give  the  more  traditional  definition. 

Let 
(1.5.1)    x^  <  X2  <  .  .  <  Xj^ 

be  a  sequence  of  strictly  increasing  real  numbers  called 
the  knots  of  the  spline  function.   We  may  say  Sj^(x)  is  a 
spline  function  of  degree  m  having  the  knots 

x-L,  ^2'  '    •     '     ^k 
if  it  satisfies 

a)  Sjj,(x)  e  C"""^  (-»,«')  ; 

b)  In  each  interval  (x^,    Xj^^-^)  ,  including  (-<^  ,x-^) 
and  (xj^,°°),  the  restriction  of  Sj^(x)  to  (Xj_,  x^_^^)    is  a 
polynomial  of  degree  at  most  m.   Thus,  a  step  function 
Sq(x)  may  be  regarded  as  a  spline  function  of  degree  0, 
while  a  spline  function  of  degree  1  is  a  polygon  (broken 
line  function)  with  possible  corners  at  some  or  all  of 
the  possible  corners  at  some  or  all  of  the  points  (1.5.1). 
Similarly,  S2(x)  has  a  graph  composed  of  a  sequence  of 
parabolas  which  join  at  the  knots  continuously  together 
with  their  slopes.   Both  for  a  smoother  approximation  and 
for  a  more  efficient  approximation,  one  has  to  go  to 
piecewise  polynomial  approximation  with  higher  order 
pieces.   The  most  popular  choice  continues  to  be  a 
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piecewise  cubic  approximating  function.  Various  kinds  of 
cubic  splines  are  in  use  in  numerical  analysis.  The  ones 
most  commonly  used  are  complete  cubic  splines,  periodic 
cubic  splines  and  natural  cubic  splines. 

A  spline  function  of  degree  m  with  k  knots  is  repre- 
sented by  a  different  polynomial  in  each  of  the  k  +  1 
intervals  into  which  the  k  knots  divide  the  real  line.  As 
each  polynomial  involves  m  +  1  parameters,  the  spline 
function  involves  a  total  of  (m  +  1)  (k  +  1)  parameters. 
However,  the  continuity  conditions  stated  earlier  impose 
certain  constraints  on  those  parameters.  At  each  knot, 
the  two  adjoining  polynomial  arcs  must  have  equal 
ordinates  and  equal  derivatives  of  order  1,2,..  ., 
m  -  1.  Thus,  m  constraints  are  imposed.  It  is  easy  to 
see  that  every  spline  function  s (x)  of  degree  m  with  the 
knots  X-,  ,  X2/  .  .  ,  X]^  has  a  unique  representation  in  the 
form 

k 

(1.5.1)  s(x)  =  Pj_(x)  +   ^  c.  (X  -  x.)^ 

j=l  ^       J 

where  Pr^(x)  denotes  a  polynomial  of  degree  m  and 

(1.5.2)  x^"^  =  x'"      X  >  0 

=   0       X  <  0  . 
Also 

(1.5.3)  Cj  =  (l/(m)!)  [  s("')(Xj  +  )  -  S^^U^y)     }  . 

The   class   of    "natural"    spline    functions    was    intro- 
duced   by    Prof.    Schoenberg     [1946].       A    spline    function    s (x) 
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of  odd  degree  2p-l  with  knots  x-^r    ^2'    '    ''    ^k  ^^    called 

a  natural  spline  function  if  the  two  polynomials  by  which 

it  is  represented  in  the  two  end  intervals  (-   ,x-j^)  and 

(Xi^,+   )  are  of  degree  p-1  or  less.   It  is  easy  to  express 

the  natural  spline  functions  by 

k  -  . 

(1.5.4)    s(x)  =  Pp-i(x)  +  I      C.  {x-Xj)^^P"^ 


D=l  ' 


where 


I      C-    x-^  =0,   r  =  p,  p+1,  .  .,  2p-l. 
j=l   ^   ^ 

The  following  theorem  states  an  important  interpola- 
tion property  of  natural  spline  functions. 

Theorem  1.9  Let  (x^,  y^^),  i=  1,  2,  .  .  ,  k,  be  given 
data  points,  where  the  x^'s  form  a  strictly  increasing 
sequence,  and  let  p  be  a  positive  integer  not  exceeding  n. 
Then  there  is  a  unique  natural  spline  function  s(x)  of 
degree  2p  -  1  with  the  knots  Xj_  such  that 
(1.5.5)    s(Xi)  =  y^  ,   i  =  1,  2,  .  .  .  ,  k. 

Natural  spline  functions  possess  certain  impressive 
optimal  properties  and  can  be  shown  to  be  the  "best" 
approximating  functions  in  a  certain  sense.  This  is  the 
content  of  the  next  theorem. 

Theorem  1.10   Let  P(x)  be  the  unique  natural  spline 
function  that  interpolates  the  data  points  {x^,Y^)t 
i  =  1,  2,  .  .  ,  k,   in  accordance  with  Theorem  1.7.    Let 
f  (x)  be  any  function  of  the  class  C^P'  that  satisfies  the 
conditions 
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(1.5.6)  f(x^)  =  y-L  ,    i  =  1,  2,  .  .  ,  k. 

Let  (a,b)  be  a  finite  interval  containing  all  the  knots 
Xj.  Then 

b  b        ^ 

(1.5.7)  /  [f  (P)  (x)]2  dx  >     /[s'PMx)]''  dx 
a  a 

with  equality  only  if  f(x)  =  s (x) . 

The  effectiveness  of  the  spline  approximation  can  be 
explained  to  a  considerable  extent  by  its  striking  conver- 
gence properties.   Interesting  contributions  were  made  by 
J.  N.  Ahlberg  and  E.  N.  Nilson  [1964],  C.  DeBoor  and  G. 
Birkhoff  [1964],  A.  Sharma  and  A.  Meir  [1967],  M.  J. 
Marsden  [1972],  T.  R.  Lucas  [1974],  E.  W.  Cheney  and  F. 
Schurer  [1968],  C.  A.  Hall  [1968],  C.  A.  Hall  and  W.  W. 
Meyer  [1976],  and  A.  K.  E.  Atkinson  [1968].   As  a  good 
reference  on  splines  which  offers  a  good  comparison  of  the 
approximating  properties  of  polynomials  and  splines,  we 
recommend  A  Practical  Guide  to  Splines  by  C.  DeBoor 
[1978] . 

First  we  discusss  error  analysis  for  the  class  of 
functions  f(x)  6  C^  '  with  period  one.   Let 

(1.5.8)  Cx.}J^o:  0  =  x^o  <  ^nl  <  •  •  •  <  ^n,kn  =  ^ 
be  a  division  of  [0,1]   of  mesh  gauge 

(1.5.9)  h^  =  ^^^0<i<]^    (hn,i)  ' 

—  —  n 
where 

h      =  V      —    ■a        ■ 

n,i    n,i    n,i-l  ' 
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A  periodic  cubic  spline  function  Y^(x)  is  .  a  function 
composed  of  a  cubic  polynomial  in  each  of  the  intervals  of 
{Xj^}^_Q  with  the  requirement  that 

y^(x)  6  C^^^  [0,1] 
and 

Yn^^^  (0)  =  Yn^"-^  (D  .  i  =  0,  1,  2. 
It  was  observed  by  Walsh,  Ahlberg  and  Nilson  [1962]  that 
there  exists  a  unique  periodic  spline  function  Y^(^)  which 
interpolates  f  (x)  at  the  points  x^^  ■^.  It  was  shown  that 
y  (x)  and  y'^(x)  converge  uniformly  to  f{x)  and  f'(x) 
respectively  as  h^^  ->  0.  Later  Ahlberg  and  Nilson  [1966] 
studied  the  more  delicate  question  of  the  convergence  of 
y"^(x)  to  f "  (X).  Writing 
(1.5.10)       A^^i    =        hj,^i^i/(h^^i   +  hn,i  +  l)    ' 

i    =    1,     2,    .    .     ,    k^    -    1 
and 

^n  =  "'^^Olilk  l^n,i  -1/2  1, 
where  for  m  =  k^,    X^   ^^-^      is  taken  as 

(hn,l  ^  hn,m)/hn,i  , 
they  show  that 

y"^(x)   ->  f"(x) 
uniformly  provided  that 

h^  ->  0  and  A^  ->  0  . 
After  this  result,  I.  J.  Schoenberg  [1964a]  raised  the 
question  that  it  would  be  very  interesting  to  find  out  to 
what  extent  the  condition  Aj^  ->  0  is  really  necessary  in 
the  above  mentioned  theorem.   The  above  theorem  together 
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with  tbe  open  problem  of  Schoenberg  lead  to  important 
contributions  by  Birkhoff  and  DeBoor  [1964] ,  and  Meir  and 
Sharma  [1969]  which  we  turn  to  describe. 

In  1964,  Garrett  Birkhoff  and  Carl  DeBoor  made  the 
following  contribution.   Let  f(x)  6  C'[0,1]  and  let 
(1.5.11)   [x^l^^o'   0  =  Xq  <  x^  <  .  .  <  Xj^  =  1 
be  a  partition.   The  function  f(x)  is  now  interpolated  by 
a  cubic  spline  function  s (x)  (called  a  complete  cubic 
interpolation  spline  function)  which  means  that  s (x)  is  a 
cubic  polynomial  when  restricted  to  each  interval 
(Xj_  ,Xj^^-]_)  ,  and  s(x)  6  C^^^[0,1].    Moreover  s(x)  is 
uniquely  defined  by  the  conditions 
(1.5.13)   f(xj_)  =  s(x^)  ,       i  =  0,  1,  .  .  ,  k  ; 

f  •  (0)  =  s'  (0)  , 

f  (1)  =  s-  (1)  . 

This  first  important  result  concerning  the  error  analysis 

yielded  the  following  theorem. 

Theorem  1.11   Let  f(x)  G  C^'*)[0,1].   Denote 
^(r)  ^  f(r)_  3(r)  . 

There  are  constants  Cj.(m),  r  =  0,  1,  2,  3,  depending 

only  on  m  >  0,  such  that 

(1.5.13)   le^^Nx)  I  <_  c^(m)  h"*"^  ||f(^)||  , 

r  =  0,  1,  2,    3, 
provided  that 
mr^  <_  m  , 

'^i  "  ^i  +  1  ~  ^i  ' 
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h  =  max  hj  , 

itij^  =  [max(hj_)  ]  /  [min(hj_)  ] 
and  II    II  denotes  the  supremum  norm. 

The  authors  go  a  step  further  and  prove  a  convergence 
theorem  related  to  f  6  C^^Ho,!]. 

Theorem  1.12   Let  f"'(x)  be  absolutely  continuous  on 
[0,1].   Let  {xj}^_Q,n  (where  k  depends  on  n)  be  a  sequence 
of  partitions  of  [0,1]  such  that  h^^  =  maXj_hj^  j^  ->  0  as 
n  ->   .   Let  mj^  ^  <_  m  as  n  ->   .   Let  e^^lx)  be  the  error 
incurred  when  f(x)  is  interpolated  by  a  spline  function  on 

CXj^}i  =  o'"*   ^^^" 

uniformly   on    [0,1]    as   n   ->  °°    . 

The  next  important  development  came  with  some 
interesting  results  by  Prof.  A.  Sharma  and  A.  Meir  [1967] 
concerning  degree  of  approximation  of  spline  interpola- 
tion. This  paper  does  away  with  some  annoying  assumptions 
under  which  uniform  convergence  of  the  interpolating  cubic 
spline  and  its  derivatives  was  proven  earlier  (see  above 
for    these    restrictions). 

Theorem    1.13       Let    f(x)    be    continuous    and    periodic    with 
period   unity.      Let 
(1.5.15)       q^^    =    max^^j     ^K,i/K,j^ 
where 

h 


n,i       ■^n,i  +  l  n,i 


Let  Sj^(x)  be  the  cubic  spline  of  period  unity  with 
joints  (or  knots)  x^  • ,   i  =  0,  1,  .  .  ,  n  in  [0,1],  such 

11  ,  X 
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that  s  (x)  interpolates  f (x)  at  the  joints.   Let 

||g||  =  maxj^|g(x)|   for  g  6  C[0,1]  , 
and 

co(g,h)  =  max  C  |g(u)  -g(v)|  :  |  u-v  |  <_  h  },  h  >_  0  . 
The  authors  prove 

i) 

(1.5.16)  ||f  -  s^ll  <  (1  +  q^^)    (o(f,h^)  ; 
ii)  if  f  e  C^^^ ,  then 

(1.5.17)  ||f(^)   -  3(^)^11   <  76  hj^l-^  a)(f',hn)  , 

r  =  0  ,  1  ; 
iii)  if  f  6  C^^^ ,  then 

(1.5.18)  [  |f  (^)  -  s(^)j^|  I  <  5  h^2-r  ^  (f..,h^)  , 

r  =  0,  1,  2  ; 
iv)   if  f  e  C<^^ ,  then 

(1.5.19)  I  |f  (^)  -  s^^^nM  <  C  h^^-^   aj(f' ",h^)  r 

r  =  0,  1,  2,  3  ; 

where 

2 


C  =  1  +  q^(l.q^) 


C  =  1  +   (1  +  P^)2/(2  -  P^) 
with 

Pj^  =  maxj^   (^n,i/^n,j^  ^°^      ^    "  ^~^'  ^"^^ 
satisfying 

P   <  2  . 

From  these  results  one  can  draw  the  obvious  conclu- 
sions regarding  uniform  convergence  of  the  interpolating 
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splines  and  derivatives.  The  arguments  are  surprisingly 
simple.  The  uniform  convergence  of  s^"  to  f",  which 
follows  from  iii) ,  had  been  proved  earlier  by  Ahlberg  and 
Nilson  (see  above)  under  the  additional  assumptions  that 
the  mesh  become  eventually  uniform,  i.e., 
(1.5.20)   lim„_>^  [h^,i/(hn,i^h^,i.l)]  =  1/2  • 

Parabolic  Spline  Interpolation 
Many  interesting  results  were  obtained  by  M.  Marsden 
[1974]  concerning  the  approximation  of  functions  by  even 
degree  splines.  Of  particular  interest  are  the  simple 
parabolic  splines.  If  break  points  are  the  same  as  the 
interpolated  points,  then  the  resulting  spline  is  ill- 
behaved,  as  can  be  seen  by  simple  examples  (DeBoor 
[1978]).  On  the  other  hand,  if  we  take  the  interpolated 
points  midway  between  break  points,  the  parabolic  splines 
are  very  well-behaved.  In  fact  in  the  first  theorem  given 
below,  a  good  approximation  to  a  continuous  function  is 
assured  with  no  conditions  on  the  partition  other  than  the 
length  of  the  largest  subinterval  being  small. 

We  first  give  some  necessary  notation.   Let 

(1.6.1)  fx^]^^^o=  0  =  Xq  <  x^  <  .  .  <  Xj^  =  1 
be  a  fixed  partition  of  [0,1].   Set 

(1.6.2)  ^i  ~  ^i  ~  ^i-1  '    ^  ~  "i^^ih.  , 
z^  =   (Xj^  +  x^_j_)  12    , 

^0  =  ^\  '   ^i  =  h.^^L/^'^i  ^  ^i  +  l)  ' 
Cj_  +  a-j_  =  1  ,   for  i  =  1,  2,  .  .  ,  n. 
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Let 

y  e  C[0,1]  ,   y(0)  =  yd)  , 
I |y| I  =  sup  C  |y(x) |  :  0  <  x  <  1  3 
such  that  y  is  extended  periodically  with  period  1. 

A  function  s (x)  is  defined  to  be  a  periodic  quadratic 
spline  interpolant  associated  with  y  and  fx£}j|_Q   if 

(1.6.3)  a)  s(x)  is  a  quadratic  expression  on  each 

(Xj^_^  /  Xj^ )   , 

b)  s(x)  e  C  [0,1]  , 

c)  s(0)   =  s(l)  ,  s' (0)  =  s' (1)  , 

d)  sCzj^)  =  y(Zj_)  ,   i  =  1,  2,  .  .  ,  n. 

The  following  theorems  were  obtained  by  Marsden. 
Theorem  1.14  (Marsden).   Let  fx^}^_Q  be  a  partition  of 
[0,1],  y(x)  be  a  continuous  1-  periodic  function  and  s  (x) 
be  the  periodic  quadratic  spline  interpolant  associated 
with  y  and  Cx^}"_q. 
Then 

(1.6.4)  ||s.||  <  2  i|y||  ,  MsM  <  2  ||y||  , 
I  l^ii  I  <  2  to(y,h/2)  , 

I  |e|  I  <  3  a3(y,h/2)  . 
(where   s^  =  s(x-)  and  e^^  =  y(Xj_)  -  s(Xj^)  ). 
The  constant  2  which  appears  in  the  first  of  the  above 
equations  can  not,  in  general,  be  decreased. 
Theorem  1.15  (Marsden).   Let  y  and  y'  be  continuous  1- 
periodic  functions.   Then 

(1.6.5)  I  |s'^|  I  <  2|  |y'  I  I  , 
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I  le'il  I  <  3  I  |y'  I  i  , 

I  le^l  I  <  h  aj(y'  ,  h/2)  , 

lleill  <  h  IIy'II  , 

I |e  II  <  (5/4)  h  I |y' II  , 

I  le'ij  I  <  3  a)(yMi/2)  , 

!|e'  I  I  <  (9/2)  (^(y',h/2)  , 

I  |e    I  I    <       (13/8)h      a3(y'  ,h/2)     . 
Theorem   1.16      (Marsden).   Let  y,    y',    and  y"   be   continuous 
1-   periodic    functions.      Then 
(1.6.6)  I  |e^|  I    <    (1/8)    h2   aj(y'  '  ,h)     , 

I  |e'^|  I    <    (1/2)    h     aj(y"  ,h)     , 

Me'll   <  2  h  llyjl   , 

I  |e|  I    <     (5/8)    h2    ||y"|  I     , 
|e"  (X)  I    <    [1    +    (h/h^)]    co(y",    h)     , 
^i    <    X   <    x.^^    . 

Theorem    1.17    (Marsden).    Let      y,    y',    y''r    and    y'''    be 

continuous    1-   periodic    functions.    Then 


(1.6.7)  I |e| I  <    (17/96)    h^    | |y 


III 


e 


<     (11/24)    h2     I  |y'  "  I  I     , 


e"  I  I         <    [h^    +    (2   h^/s   h^)]     I  |y' 


^i  +  1    • 


Optimal  Error  Bounds  for  Cubic  Spline  Interpolation 
An  interesting  application  of  the  theorem  of  Birkhoff 
and  Priver  [1967]  (discussed  above)  was  given  by  Hall 
[1968]  and  subsequently  by  Hall  and  Meyer  [1976],  concern- 
ing optimal  error  bounds  for  cubic  spline  interpolation. 


28 


(4) 


In  order  to  describe  these  results  let  f  G  C^^'  [0,1]  and 

let  s(x)  be  the  complete  cubic  spline  function  satisfying 

the  conditions  (1.5.13).   The  main  result  of  Hall  and 

Meyer  may  now  be  stated. 

Theorem  1.18  (Hall  and  Meyer).   Let  s(x)  be  the  unique 

complete  cubic  spline  interpolation  satisfying  (1.5.13). 

Suppose 

fee'"*)  [0,1]  . 
Then  for  0  _<  x  _<  1 
(1.7.1)    If (^) (X)  -  s(^>  (X) I  <  c^  h^-^  [if (4)  i  I 

r  =  0,  1,  2 
with 

h  =  max(x-^j^  -  x^)  ,    Cq    =    5/384  , 
c-|_  =  1/24  ,   C2  =   3/8  . 
Further,  the  constants  Cq  and  c-|_  are  optimal  in  the  sense 
that 

'l-'7-2)    ^r  =  sup  l|(f  -  s)^^Ml 

where    the    supremum    is    taken   over    all    Cx.}j^_q    partitioning 
[0,1]     and      over    all    f    e    c('*)[0,l]     such    that    f^^^     is    not 
identically   equal    to    zero. 

Varma    and    Katsifarakis     (in    press)     were    able    to 
resolve    the    cases    of     f     e    C^^^     and     f     e    C^^^     in     the 
following   theorems.      Let    s (x)    be   the   unique   complete   cubic 
spline    satisfying    the    relationship: 
(1.7.3)         s(x^)     =    f(x^)     ,    i    =    0,    1,     .     .     ,    k    ; 

s' (x^)     =    f ■ (x^)     ,    i    =    0,    k. 
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Theorem  1.19  If  f,  f,  f ',  and  f"'  are  continuous  on 

[0,1] ,  then 

(1.7.4)    |s(^)  (X)  -  f  (^^  (X)  I 

1  ^r   ^^"^  '"^^0<x<ll^'  '  '  (^H 
r  =  0,  1,  2 

where 

Cq  =  1/96  +  1/27  ,  c^  =  4/27  , 

C2  =  1/2  +  4/ (3  /3)  . 
Theorem  1.20   If  f,f',  and  f"  are  continuous  on  [0,1], 
then 


(1.7.5)    |s^^Mx)-f^^Mx)|<a^  h^"-^  w(f  '  '  ,h] 


i^(i^)  i^\    -    ^(^)  1^)  I    ^    ^      h2-r 
where 

aQ  =  13/48  ,  a-L  =  5/6  ,  82  =  4 


CHAPTER   TWO 

BEST         ERROR      BOUNDS      FOR      DERIVATIVES      OF 

TWO      POINT      LIDSTONE    POLYNOMIALS 


Introduction  and  Statement  of  Main  Theorem 
Let  u  e  C^'^[0,h]  be  given  and  let  V2^_j_  be  the  unique 
Hermite  interpolation  of  degree  2m  -  1  matching  u  and  its 
first  m-1  derivatives  u^^'  at  0  and  h.  Let  e  =  V2j^_]_  -  u 
be  the  error  function.  For  the  special  cases  m  =  2  and 
m  =  3,  G.  Birkhoff  and  A.  Priver  [1967]  obtained 
pointwise  optimal  error  bounds  on  the  derivatives  e^  ', 
0  <^  k  <_  2m  -  1  in  terms  of  h  and  maxQ^ ^^^ j^  |  u  ^^"'Mx)  |  . 
These  results  are  described  in  detail  in  Chapter  One. 
Birkhoff  and  Priver  note  that  for  the  cases  m  >  3,  their 
method  is  not  likely  to  give  analytically  exact  bounds, 
though  it  can  be  adapted  to  give  numerical  approximations 
to  pointwise  exact  error  bounds.  In  the  next  chapter,  we 
will  directly  apply  the  results  of  Birkhoff  and  Priver  to 
the  case  of  u  in  c  ^  ^  "^ '  [0,h]  and  the  inter  po  1  a  tory 
polynomial  W2j^_i  which  matches  u  at  0  and  h  and  which  also 
matches    the    2nd    through    mth   derivatives    of   u   at    0    and   h. 

Analogously  to  using  Hermite  interpolation 
polynomials,  one  may  choose  to  approximate  a  given 
function  u(x)  in  C^'^LOfh]  by  the  so-called  Lidstone 
interpolation   polynomial    L2^_]_[u,x]    of    degree    <_   2m    -    1 
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matching  u  and  its  first  m  -  1  derivatives  u^^^'  at  0  and 
h.   Thus  L2^_;i_[u,x]  satisfies  the  following  conditions 
(where  we  assume  h  =  1)  : 

(2.1.1)  L2m-1^^^^  [^'  °1   =u(2p)  (0)  , 

L2j,_i<2p)[u,   1]  =  u(2p)(i)  , 

p  =  0,  1,  .  .  ,  m  -  1. 
The  explicit  formula  for  L2[p_;l[u,x]  is 

(2.1.2)  L2m-lf^'^^  =   ^  ^     (^>   i^""^ 

i  =  0 

m-l  ,_., 
+  I    u(2i) (0)   ^(1-x) 
i  =  0 
where 

(2.1.3)  A.  (x)  =   2^  Bojxi  (1+x)  ,  for  i  >  1 

(2i+l)!       2 

and 

(2.1.4)  AqIx)  =  X  . 

Here  B  (x)  denotes  the  Bernoulli  polynomial 
n 

(2.1.5)  B^(x)  =  I   (J)  x^  B^_j^ 

and  where  the  constant  B-  is  given  by 

(2.1.6)  B.  =  Hi)     Bj^  ,  Bq  =  1  . 

-'         k  =  0 

That    (2.1.2)    in  fact    satisfies    (2.1.1)    follows    from 
the   facts 

a|2p) (0)    =  0    ,       p   =    0,    1,    .    .    ,    i    ; 

(2.1.7)  A|2p)  (1)     =  0    ,       p    =    0,    1,     .     .     ,    i    -    1    ; 
A|2i) (1)     =  1    . 

The     main     object     of     this     chapter     is     to     obtain 
pointwise  optimal   error  bounds   for 
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elj)  (X)  =f  (J)  (X)  -  L2n,_|^^  [f,x] 
in  terms  of  U  =  maxQ<. ^^^  ]_  |  u  ^^"'Mx)  |  .  Here  L2n,_p  ^  [  f  ,x] 
denotes  the  j th  derivative  of  the  Lidstone  polynomial 
defined  by  (2.2.2).  An  important  role  in  Theorem  2.1  (see 
below)  is  played  by  the  Euler  polynomial  Q2m(^)  °^  degree 
2m  given  by  the  formula 

1 

(2.1.8)  Q2m(x)    =    -  G^L^^'t)    Q2m-2(^^'^^    ,    m   =    1,    2,    .    . 

where 

(2.1.9)  Qo(x)    =    1 
and 

(2.1.10)  G^(x,t)     =t(x-l),       0<t<x<l 

=  x(t-l),Oj<x<t<l. 
We  may  now  state  the  main  theorem  as  follows. 
Theorem  2_^-   Let  u(x)  8  C^^'lO,!]  and  let  L2jj,_i[u,x]  = 
L2^_-l(x)  be  the  unique  polynomial  of  degree  <_  2m  -  1 
satisfying  the  conditions  (2.1.1).   Then,  for  0  <_  x  <_  1 , 
with 

U  =  maXo^^<;^|u<2m)  (x)  |  , 

(2.1.11)  |u(2J)(x)  -  L2ni_|2D)  (X)  I  <  u  Q2m-2j(^)  ' 

j  =  0,1,  .  .  ,  m  -  1 

<  U  Q2m-2j(l/2)  ' 

j  =  0,  1,  .  .  ,  m  -1 
and  for  j  =  1,  2,  .  ,  ,  m 

(2.1.12)  |u(2j-l)(x)  -  L2j,_j2J-l)  (X)  I 
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<  U  C(l-2x)  Q2m+2-2j'(^) 

<  U  lQ2m+2-2j'(0H 

where  for  a  given  integer  k,  Q2k^^^  ^^  ^^^  well  known 
Euler  polynomial  defined  by  (2.1.7).  Moreover,  (2.1.11) 
and  (2.1.12)  are  both  best  possible  in  the  sense  that 
there  exists  a  function  u(x)  e  C^"^[0,1]  such  that  (2.1.11) 
and  (2.1.12)  become  equality  for  every  x  S  [0,1]. 

From  (2.2.11)  and  (2.1.12)  follow  immediately  the 
also  exact  bounds 

(2.1.13)  l|u(2J)  -  L2^_|2J)||  <  Q2^_2j(l/2)  ||u(2m)|,  , 

j  =  0,  1,  .  .  ,  m  -  1 
and 

(2.1.14)  ||u(2J-l)  _  L2„_|2J-l)i| 

<  lQ2m.2-2j'(0)l  M^^'^Ml  . 

j  =  1,  2,  .  .  ,  m  -  1 

where  | |   | I  denotes  the  supremum  norm  on  [0,  1] . 


Preliminaries 
It  is  well  known  that  the  Bernoulli  polynomials 
defined  by  (2.1.5)  satisfy 

(2.2.1)  B^' (X)  =  nB^_^(x) 
and 

(2.2.2)  B^(l-x)  =  (-1)"  B^(x)  . 
In  particular  it  follows  that 

(2.2.3)  B2n+i(l/2)  =  0  . 
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From  (2.2.1),  (2.2.3)  and  (  2  . 1 .  3 ) -  (2  . 1 .  6)  ,  we  obtain 

(2.2.4)  A^'Mx)  =  Ai_i(x)  ",  A^(0)  =  0  ,  A^(l)  =  0  , 

i  >  1  . 
The  proof  of  Theorem  2.1  depends  on  repeated  use  of 
the  kernel  G^(x,t)  defined  by  (2.1.10).   Let  us  consider 

1 

(2.2.5)  g(x)  =  /  G,(x,t)  r(t)dt 

0 
X  1 

=  /  (x-l)t  r(t)dt  +  /  (t-l)x  r(t)dt  . 
0  x 

On  differentiating,  we  have 

X 

g*  (x)  =  /  t  r(t)dt  +  (x-l)x  r(x) 
0 

1 
-   /  (t-1)  r(t)dt  -  x(x-l)  r(x) 

X 

x  1 

=  /  t  r(t)dt   +  /  (t-1)  r(t)dt  . 
0  X 

Differentiating  once  more  with  respect  to  x  we  obtain 

(2.2.6)  g''(x)  =  X  r(x)  -  (x-1)  r(x)  -  r(x)  . 
Also 

(2.2.7)  g(0)  =  g(l)  =  0  . 

Let  r(t)  =   -/^ni-l^^^  ^^     (2.2.5).   From  the  above 

discussion  it  follows  that 

1 
g(x)  =  /  G^L^X't)  Aj^_^(t)dt 
0 

satisfies 

(2.2.8)  g"  (X)  =  A^_^(x)  ,  g(0)  =  g(l)  =  0  . 
From  (2.2.4)  we  also  know  that  for  i  >^  1 

Ai"(x)  =   A^_^(x)  ,   A^(0)  =  0  ,   A^(l)  =  0  . 


35 


Therefore 

1 

(2.2.9)  g(x)  =   A  (X)  =   /g^Ix,!)  Aj^_^(t)dt.. 

0 

From  (2.1.9)  it  follows  that 

(2.1.10)  G^(x,t)  <_  0  . 

Also       AQ(t)     =t>_0,     0<_t<_l.       Therefore    we 
obtain    from     (2.2.9)     that 

(2.2.11)  A;l(x)    10,      0    <_  x  £   1    . 

On  using  (2.2.9),  (2.2.10),  and  (2.2.11),  we  can  assert 
that 

(2.2.12)  A2(x)  >    0    ,    0  _<  X  £  1  . 

Inductively,  it  follows  that  A^(x)  >_  0  for  0  <_  x  <_  1 
provided  m  is  an  even  positive  integer  and  Aj^(x)  £0  , 
0  <_  X  <_  1  if  m  is  an  odd  positive  integer.   This  property 
of  A  (x)  will  be  needed  many  times  in  the  proof  of  the 
theorem. 

The  following  iteratively  defined  kernels  comprise 
the  essential  machinery  of  the  proof.   Define 

1 

(2.2.13)  G2(x,t)  =    /  Gi(x,y)  Gi(y,t)dt 

0 

and  inductively 

1 

(2.2.14)  G^(x.t)  =    /  G-^(^.Y)    G^_j^(y,t)dy   n  =  2,  3,  .  . 

0 

From    (2.2.10)    and    (2.2.13)    it  follows  that 

(2.2.15)  G2(x,t)     >_   0    ,       G3(x,t)     <    0    ; 

0<x<l,       0<t<l. 
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In  general 

(2.2.16)  (-l)^Gj^(x,t)  >  0  ; 

Oj<x_<l,   0£t_<l. 
Finally,  let  us  define 

1 

(2.2.17)  h(x)  =    /  Gj^(x,t)  q(t)dt  . 

0 

We  note  again  that  h(x)  uniquely  satisfies 

j^(2n)  (X)  =  q(x) 

(2-2.18)    ,„,.        ,.,. 

h(2k) (0)  =  h(2k)(i)  =  0  ,   k  =  0,  1,  .  .  ,  n-1  . 

We  also  need  some  of  the  known  properties  of  Euler 

polynomials  introduced  in  (2.1.7)  and  (2.1.8).   We  can 

easily  verify  that 

Q2A'<^)  =  Q2n-2(^) 
(2.1.19) 

Q2n(0)  =  Q2n(l)  =  0  • 
Furthermore, 

Q2ii^PM0)  =  02^2?)  (1)  =  0  ,  p  =  0,  1,  .  .  ,  n-1, 
(2.2.20)   Q2ji2^Ml)  =  Q2il^'^N0)  =  (-1)" 

Q2ji2J)(x)  =  (-1)3  Q2n-2j(^) 

Using  (2.2.13)  we  note  that 

1 
Q2(x)  =  -    /  Gj_(x,t)  dt  , 
0 

1 
Q^(x)    =    -      I    Gi (x,t)  Q2(t)  dt 
0 

1  1 

/  G.(x,t)  [    /  Gi(t,y)  dy  ]  dt 
0  0 

1 
/  G2(x,t)dt  , 


and  in  general. 
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1 
m 


(2.2.21)   Q2m(x)  =  ("D"'    /   G^(^'t)  dt   . 

Explicitly  some  of  the  first  Euler  polynomials  are 

given  by 

Qp(x)  =  x(l-x)  ,  Q.(x)    =   x^(l-x)^+x(l-x)  , 
2l  41 

Qg(x)  =  x^(l-x)^+3x^(l-x)^+3x(l-x)  . 
6! 

Proof  of  Theorem  2.1 
Let  P^m-l  '^s'^ote  the  class  of  polynomials  of  degree 
<_  2m-l.   Following  the  notation  used  by  Birkhoff  and 
Priver  [1967]  we  shall  denote 

9i"':iG„(x,t) 


(2.3.1)    G^^i'^i)  (x,t)  = 


'3x^  3t^ 

Since   L2j^_]^[u,x]    =   u(x)    for  u(x)    G    P2ni-1   ^^   follows 
from    the   Peano   theorem   that   for  u  e   C^™[0,1] 

(2.3.2)  e(x)     =:    u(x)    -    L2m-lf"'^^ 

=      J    G^(x,t)    u<2m) (t)    dt 
0 

where  G  (x,t)  is  the  Peano  kernel  defined  by  (2.1.10)  and 

(2.2.14).   Differentiating  (2.3.2)  we  have 

(2.3.3)  e(2j)(x)  =  u(2J)(x)-  L2n,_  j^^  ^  [u,x] 

=    /  G^(2j,0)(^^t)  u(2ni)(t)  dt  . 
0 

Let    us     substitute    u(x)     =    Q2m^^^     ^^^    defined    by 
(2.1.7))    in    (2.3.3)    and  use  various    properties    as   given  by 
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(2.2.20)  and  (2.2.21).   We  then  obtain 

(2.3.4)  Q2rJ^^^(^)  -  L2n,-l^^^[Q2m'^l 

=   /^G^2j,0)(j,^t)Q2^2m)(t,dt  . 

We  know  from  (2.2.20) 

(2.3.5)  Q2i,^"'Nt)  =  Q2i{,^P^0)  =  (-1)""  • 
Moreover, 

(2.3.6)  Q2ii^^^(0)  =  Q2m^^^  (1^  =  0 ,  p  =  0 ,  1 ,  .  .  ,  m-1  . 
It  follows  that 

(2.3.7)  L2„_i[Q2^,x]     =    0 

identically.       Thus     (2.3.4)    can    be    rewritten    as 

(2.3.8)  Q2m^^^(x)     =    (-l)""  /     GjJ^J/O)  (x,t)    dt    . 

Next  we  note  from  (2.2.14)  that 

Gii''°^  =  G^-l(^'t)  . 
Hence 

G^4,0)(^^t)     =    G„_|2.0)(^^t)     =    G^_2(x,t) 
and   in   general, 

(2.3.9)  GjJ2J/0)  (j^^t)     =Gj^_j(x,t)     . 
From    (2.2.16)    and    (2.3.9)    we   have 

(2.3.10)  (-1)"^"^    Gj{,2J'0)  (x,t)     =     (-1)^"^    Gj^_j(x,t)     >    0 
in   the   unit   square    0_<x_<l    ,    0_<t_<l. 

Combining     (2.3.3),      (2.3.9),      (2.3.10),      (2.2.19),     and 
(2.3.8),    it  follows  that 

le(2J) (X) 1    <    U         /    |G^2j,0) (^^t) |    dt 
0 


39 


=  U  I  /  GjJ^J'O)  (x,t)  dt  I 
0 

=  "  Q2m-2j(^)  ' 

j=0,  l,..,m-l. 

This  proves  (2.1.10). 

We  next  turn  to  prove  (2.1.11).   Due  to  (2.3.9),  it 

is  enough  to  prove  (2.1.11)  for  j  =  1.   From  (2.2.14),  it 

follows  that 

X 


2.3.11)   GjJ^'O)  (x,t)  =     /  y  Gj^_T  (y,t)  dy 

0 


1 
/  (y  -  1)  G^_i(y,t)  dy  . 

X 


Therefore 

(2.3.12)     /  |g(1'°^  (x,t) I  dt 


0    " 


1    X 

<       /   /  Y  lG„_i(y,t) I  dy  dt 
0    0 

1    1 
+   /    /  (1-y)  \G^_-^(Y,t)   I   dy  dt  . 
0    X 

Recalling  (2.2.21) 

1 

Q2m-2(y)  =  (-1)"'"^   /  G^_i(y,t)dt  , 

0 

m  =  2,  3,  .  .  . 

and  the  fact  that  in  the  unit  square  0<_x£l,  0<_t<_l, 

(-I)'"-!  G^_i(y,t)  >  0  , 
we  can  assert  that 
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1 

(2.3.13)  Q2ni-2(y^  =     ^     \^m-l^^'^^\    ^^  ' 

On  changing  the  order  of  integration  in  (2.3.12)  and 
making  use  of  (2.3.13),  we  obtain 

1  X 

(2.3.14)  /  |Gjj,l'°Mx,t)  I  dt  <    /  y  Q2m-2<y)  ^y 
0  0 

1 

+  /  (1-Y)  Q2m-2(y)  ^y 

X 

=  =   X2in-2(^)  • 
Using  (2.2.20)  we  note  that 

X 

(2.3.15)  X2n,-2(^)  ="    ^  ¥02^^^^^^ 

1 

/  (1-  Y)Q2m  ^y^^y  • 

X 

On  integrating  by  parts,  we  have 

X  X 

(2.3.16)  X2ni_2(x)  =    '   Y   Q2m'  ^y^  \       ^     /  Q2m'  ^^^  ^^ 

1       1 

X       X 

=  -X  Q2m'(x)  +  (1-x)  Q2m'(><)  ^  2Q2j,(x) 

=  (1  -  2x)  Q2m'(^)  ^  2Q2^(x)  . 
Also 

(2.3.17)  X2ni-2'<^)  =  (l-2x)  Q2m'(x)  . 

Since  Q2m-2  vanishes  only  at  x  =  0  and  x  =  1,  it  follows 
that  the  critical  point  at  X2m-2(^^  inside  C 0 , 1 ]  is  only 
at  x  =  1/2.   Also  we  note  that  X2m-2(^^  "  X2m-2^°^- 
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Further 

(2.3.18)  X2^_2(l)  -  X2„_2(l/2) 

1 
/  (2x-l)  Q2m-2(y)  dy  >  0  . 
1/2 

Thus  we  conclude  that  X2m-2^^^  ^^^    ^"  absolute  maximum  at 

X  =  0  and  x  =  1.   Therefore,  from  (2.3.2),  (2.3.14),  and 

(2.3.11),  it  follows  that 

(2.3.19)  |e'(x)|  <   U    /  |Gy'°Mx,t)|  dt 

0 

1  "  X2m-2(^) 

<  U  I  (l-2x)  Q2n,'  (X)  +  2Q2j„(x)  | 

<  U  ^2m-2(l)     • 
On   using    (2.3.15)    it   follows    that 

|e'(x)|    <    U   X2^_2(l)     =    -    U   Q2„'(l)     =    U   Q2^M0)     , 

which  proves    (2.1.12). 

That  (2.1.11)  and  (2.1.12)  are  best  possible  follows 
from  the  Peano  theorem,  or  more  simply,  by  choosing  u(x)  = 
Q2p^(x),  the  Euler  polynomial  defined  by  (2.1.7).  In  view 
of  (2.2.20),  we  have  U  =:  raaxQ^^^^^  |  u'^"*)  (x)  |  =  1.  Further 
use  of  (2.2.20)  and  the  definition  of  L2,^_]_[u,x]  show  that 
L2fp_2_  [Q2mAx]  is  identically  zero.  Our  choice  of  u(x)  then 
gives  pointwise  equality  in  (2.1.11).  Similarly  it  can  be 
shown  that  (2.1.12)  is  also  pointwise  best  possible.  This 
proves    the    theorem. 

It  is  perhaps  worth  remarking  that  any  exact 
evaluation  of  the  integral  of  the  absolute  value  of  a 
Peano    kernel    results    in    an    exact    error    bound    (see    Sard 
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[1963]  or  Stroud  [1974]).  Generally  error  bounds 
resulting  from  integration  of  a  Peano  kernel  under  the 
assumption  that  u(x)  G  C^[a,b]  also  hold  for  u  having 
piecewise  continuous  kth  derivative  on  [a,b] ,  and  even  for 
u  having  (k-l)st  derivative  absolutely  continuous  on 
[a,b].  In  the  case  given  here  we  can  thus  expand  the 
class  of  functions  for  which  the  error  bounds  of  Theorem 
2.1  hold  and  hence  are  best  possible. 

As  Theorem  2.1  is  stated  for  function  u(x)  2m  times 
continuously  dif f erentiable,  it  also  holds  when  the  2mth 
derivative  is  merely  piecewise  continuous  on  [0,1]. 
Moreover  the  theorem  holds  even  for  the  case  that  u(x)  has 
its  (2m-l)st  derivative  absolutely  continuous.  In  this 
last  case  U,  instead  of  being  the  max  of  the  2mth  deriva- 
tive on  [0,1],  becomes  the  "L  infinity"  norm  of  the  gener- 
alized 2mth  derivative.  In  the  following  chapters  the 
classes  of  functions  k  times  continuously  dif f erentiable, 
the  class  of  functions  having  piecewise  continuous  kth 
derivative  and  the  class  having  k-lst  derivative  absolute- 
ly continuous  may  be  treated  as  being  interchangeable. 


CHAPTER   THREE 
MORE    POLYNOMIAL   ERROR   BOUNDS 

Introduction   and   Statement   of   Theorems 

Let    u    6    C  ^^"""^^^  [0,h]     be    given.     It    follows     from    a 

result    of    Schoenberg     [1966]    that    there    exists    a    unique 

polynomial  W2jj,  +  ]_  [u,x]    of   degree    <^   2m  +  l    satisfying 

(3.1.1)  W2n,+  i[u,0]     =    u(0)     ,  ^2m+lf^'^l    =    ^^^^     ' 
^2m+l^^^  [u,0]    =   u^P) (0)     , 

"2m+l^^^  [u,h]    =   u<P) [u,h]     , 

p=2,3,    ...,m+l. 
Theorems     3.1     and     3.2    will     give    bounds    on    u^^Mx)     - 
w^^'2rn  +  i(x)    for   the   cases   m   =   2   and   m   =   3   of   polynomials 
^2m+l    satisfying    (3.1.1). 

The  polynomial  w^  j^^  ■]_  [u  ,x]  can  be  expressed  in 
relation  to  the  Hermite  polynomial  ^2 jp_j^  [u"  ,x]  .  To 
illustrate  the  relation  between  W2j^,  +  ]_  and  V2j^_3^f  let  h  =  1 
and  let  V2~_i  [g^x]  be  the  Hermite  polynomial  of  degree  at 
most  2m  -  1  matching  g  =:  u''  and  its  first  m  -  1 
derivatives   at    0    and    1.      We   can   represent  V2j^_]_[gfX]    as 

(3.1.2)  V2n,_i[g,x]    =         AQ(x)g(0)       +  BQ(x)g(l) 

+    A^(x)g' (0)     +    B^(x)g' (1) 
+   A2(x)g' ' (0)     +    B2(x)g' ' (1) 

+    A^_-^(x)q^^-^U0}     +    B^_^(x)g("^-1)  (1) 
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where  A^(x)  and  Bj^(x),  i  =  0,   1,   .  .  ,  m  -  1  are 
polynomials  of  degree  2m  -  1  or  less  satisfying 
(3.1.3)    A^^:*)  (0)  =  6^j  ,  A^<^)  (1)  =  0  , 

j  =  0,  1,  .  .  ,  m  -  1 
B^^^)  (0)  =  0  ,  B^^^i)  (1)  =  6ij  , 

j  =  0,  1,  .  .  ,  m  -  1  . 
Define  for  i=0,  1,  ..,m-l 


1 
(3.1.4)    C^(x)    =   /  G-L(x,t)A^(t)dt  , 
0 

1 
D^(x)  =   /  Gj^(x,t)B^(t)dt  , 
0 

From  (3.1.4),  (3.1.3)  and  (2.2.5)  -  (2.2.8)  ,  it  follows 

that  for  i  =  0,  1,  .  .,  m  -  1 

(j)  rn^  =  (S.  .  .  -.  .  r.  (j) 


(3.1.5)  C^^JMO)  =  5i(j-2)  '   C-^JMl)  =  0  , 
D^(^) (0)  =  0  ,   D^(^) (1)  =  6i(j_2)  ; 

j  =  2,  3,  .  .  ,  m  +  1 
where 

(3.1.6)  C^(0)  =  D^(0)  =  C^d)  =  D^(l)  =  0  , 

and  Ci,    D^   are  polynomials  of  degree  2m  -  1  or  less. 

For  a  given  u  6  c'^"'Mo,l]  we  can  use  (3.1.5)  and 

(3.1.6)  to  give  W2j^  +  ;j_  [u,x]  in  the  form 

(3.1.7)  W2n,^.i[u,x]  =  u(0)  (1-x)   +   u(l)  X 

+  u- ' (O)Co(x)   +   u' ' (1)Dq(x) 

+  u<-^Mo)C]^  (X)  +  u^-^'  (1)D;l(^) 

+ 

+  u('^^l)  (0)Cj^_i(x)  +  u^'^*^Ul)D^_^{x) 


m 
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For  m  =  2  and  m  =  3 ,  we  give  (3.1.7)  explicitly.   For 
=  2,  if  u  e  c'°^  [0,1]  ,  then  the  unique  quintic  w^  [u,x] 


matching  u  and  its  second  and  third  derivatives  at  0  and  1 
is  given  by 

(3.1.8)  W5[u,x]  =  (1-x)  u(0)  +  X  u(l) 

+  u''(0)   [-7X/20  +  x2/2  -  x^/A    +  x^/10] 
+  u''(l)   [-3X/20  +  x^/A    -   x^/10] 
+  u'"(0)  [-X/20   +  x^/6  -  x^/6  +  x^/20] 
+  u'''(l)  [x/30   -  x^/12  +  x^/20]  . 
For  u  6  C^^^[0,1],  the  unique  polynomial  w-7[u,x]  of 

degree  <_  7,  matching  u  and  its  second,  third  and  fourth 

derivatives  at  0  and  1  is  given  by 

(3.1.9)  Wy[u,x]  =  (1-x)  u(0)  +  X  u(l) 

+    u"(0)        [-5x714    +    x^/2    -   x^/2    +    x^/2    -   x^/7] 
+    u"(l)        [-x/7    +    x^/2    -   x^/2    +    x''/7] 
+    u^^Mo)     [-13x7210    +    x^/6    -    3x^/10 
+    4x^/15    -x^/14] 
+    u^^Ml)     [4X/105    -   x^/5    +    7x^/30    -   x'^/14] 
+   u<^^  (0)     [-X/210    +    x^/24    -    3x^/40 
+    x^/20    -    x'^/84] 
+    u^"^)  (1)     [-X/280    +    x^/40    -    x^/30    +    x'^/84]     . 
The     following     theorem     concerns     the     quintic 
interpolant  w^ . 
Theorem    3.1        Let   u   6   C^[0,1]    and    let   W5[u,x]    satisfy 

(3.1.10)  W5<P) [u,0]     =    u^P^  (0)     , 

Wc^P^  [u,l]     =    u^P) (1)     ,      p    =    0,    2,    3    . 
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Denote 

(3.1.11)  e(x)  =  u(x)  -  W5[u,x] 
and 

(3.1.12)  U  =  raaxQ<j^<^|u^^Mx)  I  . 

Then  for  0<x<l,  p  =  0,  1,  2,  3,  4,  5,  the  following 
pointwise  bounds  hold: 

(3.1.13)  |e<P^  (x) I  <  U  fo,p(x) 
where 

fQ  q(x)  =  [  x^d-x)^  +  x2(l-x)2/2  +  x(l-x)/2  ]  /  6!  , 
fg^lCx)  =  [  1/60  -  x3(l-x)^/3  ]  /  4!  , 
fo,2(^)  =  [  x2(l-x)2  ]  /  4!  , 

fo,3(x)  =  [  x(x-l) (2X-1)  ]  /  12  ,    0  <  X  <  1/3 
=  [  16x^  -  105x^  +  197x  -  162 

+  66/x  -  13/x2  +  1/x^  ]  /  96  , 

1/3  <  X  £  1/2 
fQ  ^(x)  =  [  48x^  +  42x'^  -  lOOx^ 

+  54x2  _  i2x  +  1  ]  /  12(l-x)^  , 

0  £x  <  1/3 
=  [  -6(x-l/2)2  +  1/2  ]  /  12  ,  1/3  <  X  <  2/3 
fQ^^(x)  =  -(x-1/2)'^  +  3(x-l/2)2/2  +  3/16  ,  0  <  x  <  1 
and  where  fg  2  ^^'^   ^o  3  ^^®  extended  to  the  whole  of  [0,1] 
by  even  symmetry  about  1/2. 

Furthermore,  the  functions  fg^p'  P  =  0,  2,  3,  4,  and 
5  are  pointwise  best  possible.  The  functions  fo,2'  ^0,3' 
fg  4  and  fQ  5  are  those  of  Birkhoff  and  Priver  [1967]  for 
two  point  cubic  interpolation. 

That  these  functions  also  serve  as  error  bounds  in 
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the  present  case  is  a  consequence  of  the  fact  that 
Wc"[u,x]  is  the  unique  cubic  matching  u"  and  u"'  at  0 
and  1.  In  other  words  w^"  [u,x]  is  the  Hermite  cubic 
interpolation  V3[g,x]  where  g  =  u".  The  error  bounds 
given  by  Birkhoff  and  Priver  in  terms  of  ^^^q^xKI^^  Mx)  | 
are  now  expressed  in  terms  of  U  =  tn^Xg^ ^^^  ]_  |  u  ^  '  (x)  |  (as 
g^'^)    is    in   fact      u<^h. 

Denoting 
(3.1.14)       Cp   =   maxo<x<i|fo,p(^H         P   =    0,    1,    .    .    ,    5 


we   have 


11      1    ,  c-i    =    1      1    , 


*-0    ~   i^  —=-    '  1    ~   — 

6^    6!  2    6!  2'*    4! 

C-o    =     1/3    ,  C4     =       1    ,  Cr    =    1^    . 

9"  4 !  12"  2 

From    (3.1.14)    and    (3.1.13)    it  follows   that  for  every 

u  e   C^^) [0,1] 

(3.1.15)       maXQ<j^<ile(P)  (x)  |    <    Cp   U    ;    p    =    0 ,    1 ,     .     .     ,    5    . 

Remark  3.1   Note  that 

c„  =  maxn/^/i|fn  n^^Mx)!    p  =  0,  1,  .  .  ,  5 


;p  =  ""^^o<x<ll^O,0 
u(x)  =  fQ^Q(yi)    t 

(X)   =  fo,o'^)  "  ^7^^0,0'^^ 


If  we  set  u(x)  =  fg  q (^^  then  we  have 


(6)  1  =  1 


=  fo,o(^)  ^"^  "  =  '"^^0<x£ll^O,0 
By  Remark  3.1  we  see  that  for  u(x)  =  fo,o^^'  equality  is 
attained  in  (3.1.15)  for  p  =  0,  1,  2,  3,  4,  5.  The 
constants  c   are  thus  the  smallest  possible- 

The  next  theorem  gives  error  bounds  for  w-y,  analogous 
to  the  error  bounds  for  W5  given  in  Theorem  3.1. 
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Theorem    3.2    Let    u    e    c(^)[0,l],    and    let    w-7[u,x]     be    a 
polynomial   of   degree    7   or    less    satisfying 

(3.1.16)  Wy^P)  [u,0]    =   u^P)  (0)     , 

Wy^P) [u,l]    =   u^P) (1)     ,      p   =    0,    2,3,4. 
Denote 

(3.1.17)  e(x)    =   u(x)    -  v-j[u,x] 
and   denote 

(3.1.18)  U    =   maxQ<j^<j_|u(^Mx)  I     . 

Then,     for    0    <    x    <_    1    and    0    <_    p    <_    7 ,     the    following 
pointwise   bounds   hold: 

(3.1.19)  le'P) (x) I    £   U    fi,p(x) 
where 

f^^Q(x)    =    [    x'^(l-x)'*    +    (2/5)x^(l-x)3 

+   x2(l-x)2/5    +   x(l-x)/5    ]    /    8!     , 
fl^^(x)    =     (1/5) (1/8!)    -    (1/4) (l/6!)x^(l-x)^    , 
^1,2^^)     =    x^(l-x)^/6!     , 

f^    3(x)     =    x2(x-l)2(l-2x)/240    ,    0    <    x    <    2/5 
=   x2(x-l)2(l-2x)/240 

+    T^ (x-1) ^     [    lOT^x^ 

+    2T(-15x2+2x+l)     +    5x(5x-2)     ]     /    120    , 

2/5    <    X  ^   1/2 
where 
T    =    [     (3x-l) (5x+l)     +     (x-1) (-15x2+6x+l)^/2    ]     /    12 
f^    4(x)     =    x(l-x)  (5x2-5x+l)/120,       0<    x    <     (4-/'6)/10 
=    x(l-x) (5x2-5x+l) /120 


^4     [    2T;,^2,',^3_T^2. 


+    Ti^     [    2Ti'^  (2x-^-3x'^+x) 
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+  12T-l(-5x^  +  8x^-3x)/5 
+  (10x^-18x^  +  9x-l)  ]  /  12  , 
for  (4-/6) /lO  <x<    (3-/3) /6 
where 
T^  =  [  15x2  -  9x  -  (x-1)  (3x(4-5x))^/2]  /6x(2x-l)  , 
f^^4(x)  =  x(x-l) (5x2-5x+l)/120 

+  W^x  [  10w2(2x2-3x+l) 

+  4W(15x2-21x+6) 
+  5(10x2-12x+3) 


+  5(10x2-12x+3)  ]/60  , 
for  (3-/3)/6  <  X  <  (6-/6) /lO 


and  where 


W  =  [  3(l-x)(5x-2)  +  x(3(l-x) (5x-l) )^/2 
6(x-l) (2x-l) 

fi  4(x)  =   x(x-l) (5x2-5x+l)/120  , 

(6-/6) /lO  <  X  <  1/2 


■1,4 


f^  5(x)  =  (2x-l) (lOx^-lOx+l) /120 

+  Wj^^  [  20W]^2(6x2-6x  +  l) 
+  24W-l(15x2-14x  +  2) 
+  30(10x2-8x+l)  ]  /  120  , 

0  <  X  <  (4-/6) /lO 
where 

W^  =  [  -  ISx^  +  14x  -  2  -  x(3x(4-5x) )^^2  ] 
12x2  _  ^2x  +  2 
=  (2x-l) (10x2-10x+l)/120  , 

(4-/'6)/10  <  X  <  (6-/'6)/10 
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=  (2x-l) (lOx^-lOx+1) /120 

-  12^^  [  20T2^(6x2-6x+l) 

+  24T2(-15x2+l6x-3) 

+  30(10x2-12x+3)  ]  /  120  , 

(6-  6) /lO  <  x  £  1/2 

where 

T2  =  [  15x^  -  16x  +  3  -  (x-1) (-15x^+18x-3)^/^  ] 
12x2  _  ^2x  +  2 
fi  (-{x)    =  [  -  15x2  +  5x  -  1  ]  /  10 


1,6 


-  W2'^  [  W22(x-l/2) 


+  W2(15x-7)/5  +  5x/2  -  1  ]  , 
0  <  X  <  2/5 
where 
W2  =  [  -  15x  +  7  -  (-15x+6x+l)^/2  J  /  (12X-6) 

1,6'"" 


fi  f-(x)  =  -(x-l/2)2/2  +  1/40  ,  2/5  <  X  <  1/2 


f^^-7(x)  =  2(x-l/2)6  -  5(x-l/2)'^/2 

+  15(x-l/2)2/8  +  5/32  , 
0  <  X  £  1 
and  where  f^    3,  f-j^  4,  f-^  5,  and  f-j^  g  are  extended  to 
[1/2,1]  by  symmetry  about  x  =  1/2.   Furthermore,  each  of 
the  functions  fi  ^   where  p  =  0 ,  2 ,  3 ,  .  ,  7  is  pointwise 
exact. 

Setting  dp  =  maxg^^^^^  |  f -^^pix)  |  ,  we  have 

(^•^•20)  ^0  =  (iMo)  i'  ^1  =  (|)  i 

do  =  i  1  ,  d-3  =  ^ 

2    g,  26  3    30,000 
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I     _   1  A  =      -1- 

4  "  1920  '  5    120  ' 


From  (3.1.19)  and  (3.1.20)  it  follows  that  for 
0  <_  p  £  7 

(3.1.21)  niaxQ<jj<;Ll^^^^  (x)  |  <  dp  U  . 

Remark   3.2      Analogously   to   Remark    3.1,    note    that 

(3.1.22)  maxo<x<i|fi,p(x)  I    =   niaXo<x<i  I  f  i ,  o  ^^^  (^H     • 

On    setting    u     =     f^^    q(x)     it     follows     from     (3.1.22)     that 
(3.1.21)    is  exact  for    each    p. 

The    following   would    seem    to    a    natural   generalization 
of  the  Theorems    3.1   and   3.2. 

Conjecture    3.3      Let   u   G    c^^"''^^^  [0,1]    and   let  W2m  +  i[u,x]    be 
the  polynomial  of  degree  at  most  2m   +    1   matching      u      and 
its     2nd,     3rd,     .     .     ,     (m+l)st    derivatives    at    0    and    1. 
Denote 

(3.1.23)  e(x)    =   u(x)    -   W2i„+i[u,x] 
and 

(3.1.24)  U    =   maxo<j,<i|u(2m+2)  (X)  I     ^ 

Then  for  p  =  0 ,  1 ,  2 ,  we  have 

(3.1.25)  |e'P) (X) I  <  U  fm-i,p(x) 

where 

m      .  . 

fm-1  O(^)  =  f    (-I)M;)    [x"^^2^^-x] 3  /   (2m)!  ' 

"  ^'^      4  =  0         [(m+i  +  2)  (m+i  +  1)]  '    '       ^      > 

fm-l,l<^)  =  [   ^  ,   ,  -  ^ '^^ 3  /  (2m)  1  ' 

(2m+2)(2n'^M       m+1 
m 
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fm-l,2<^^  =  [  x'^d-x)'^  }  /  (2m)  !  . 
Furthermore  (3.1.25)  is  pointwise  exact,  p  =  0  and  2. 

Analogously  to  Remarks  3.1  and  3.2,  it  may  be  that 
for  every  u  e  c^^"^^^^  [0,1]  and  p  =  0,  1,  .  .  ,  2m  +  l 
(3.1.26)   maxo<^<i|e(P)  (X)  I  <  U  maXo<x<i  I  fm-i ,  0  ^^^  ^^H  • 
If  Equation  (3.1.26)  holds  then  it  is  best  possible  as  can 
be  verified  by  choosing  u  =  fj^_]_  g  ^^^  noting  that  then 
e(x)  is  the  same  as  fj^.i  q^^^*  ^°^   p  =  0,  1,  2, 

'^^^0<x<llVl,p<^n  =  "'3^0<x<ll^m-l,0^^^  (^H  • 
Hence  if  (3.1.25)  holds  then  (3.1.26)  is  true  for 

p  =  0 ,  1 ,  2.   As 

^m-l,2*^^(^)  =  [x"'(l-x)'^]/(2m)!  , 
the  conjecture  of  (3.1.26)  is  related  to  the  following 
conjecture. 

Conjecture  3.4  Let  u  e  c^^"'^  [0,1]  and  let  V2rn_i  be  the 
Hermite  polynomial  of  degree  at  most  2m-l  matching  u  and 
its  first  m-1  derivatives  at  0  and  1.   Denote 

U  =  maxQ<jj<^|u(2"^)  (x)  | 
and 

e(x)  =  V2j^_j^[u,x]  -  u(x)  . 
Then 

maxo<x<l 1^  ^^^  ^^^  I 

<  U  maxo<x<i|dP  [x"^  (1-x)  "^Z  (2m)  !  ]  |  , 
—  —  dxP 

p  =  0,  1,  2,  .  .  ,  2m-l  . 

The  results  of  Birkhoff  and  Priver  demonstrate  Conjecture 

3.4  for  the  cases  m  =  2  and  m  =  3.   Recent  work  of  Bojanov 

and  Varma  indicates  that  Conjecture  3.4  is  in  fact  true. 
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The  next  theorem  will  concern  an  interpolatory 
polynomial  which  enjoys  a  similar  property  to  that  of 
the  above  conjectures.  Let  u  e  c'  '  [0,1].  Define 
k3[u,x]     by 

(3.1.27)  k3[u,x]    =    u(0)     (l-x)(l-2x)2    +    u(l/2)     4x(l-x) 

+  u(l)  x(l-2x)2  +  u'  (1/2)  2x(l-x)  (2x-l)  . 
Then  ko [u,x]  is  the  unique  polynomial  of  degree  3  or  less 
satisfying 

(3.1.28)  k3[u,x]    =   u(0)     ,  k3[u,l]    =   u(l)     , 
k3[u,l/2]     =    u(l/2)     ,  k3'[u,l/2]    =    u'(l/2)     . 

Theorem    3.3.      Let   u   €   C^^^  [0,1] .      Denote 

e (x)    =   k3 [u,x]    -    u(x)    , 

U  =  maxQ<j^<^|u(^^x)  |  . 
Then  for  p  =  0,  1,  2,  3,  we  have 

(3.1.29)  le^P) (x) I  <  ap  U 
where 

aQ  =  1  /  (2^  4!)  ,     a;L  =  1  /  ^^^  41)  , 
32  =  (5/2)  (1/4!)  ,     a3  =  1/2  . 
That  the  a   are  the  best  possible  can  be  verified  by 

ir 


choosing 


u(x)  =  [  x(l-x) (l-2x)2  ]  /  (2^  4!) 


Due  to  the  similarity  between  the  proof  of  Theorem 
3.3  and  several  other  proofs  in  the  following  chapters,  it 
would  be  redundant  to  prove  it  here. 
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Proof  of  Theorem  3.1 
Let  u  e  C^^^  [0,1]  .   Then 

(3.2.1)  W5[u,x]  =  (1-x)  u(0)  +  X  u(l) 

+  u''(0)   [-7x720  +  x2/2  -  x'^/4  +  x^/10] 
+  u''(l)   [-3X/20  +  x'^/4  -  x^/10] 
+  u'''(0)  [-X/20   +  x^/6  -  x^/6  +  x^/20] 
+  u'''(l)  [x/30   -  x^/12  +  x5/20] 
is  the  only  polynomial  of  degree  <_  5  satisfying 

(3.2.2)  w^^P) [u,0]    =   u(P) (0)     , 

W5(P) [u,l]     =    u(P) (1)     ,  p    =    0,    2,    3    . 

Define 

(3.2.3)  e(x)    =   u(x)    -   W5[u,x]     . 
Then 

(3.2.4)  e^P^O)     =    0,    e^P^l)     =    0,    p    =    0,    2,    3, 
and 

(3.2.5)  e'^Mx)    =   Q(x)    =:    u^^Mx)     . 

In  other  words,  e(x)  is  the  unique  solution  of  the 
differential  equation  (3.2.5)  with  boundary  conditions 
(3.2.4).   We  can    rephrase    (3.2.4)    and    (3.2.5)    as 

(3.2.6)  d£e    =   y(x)     , 
dx^ 

e(0)    =    0,    e(l)    =    0    , 

and 

(3.2.7)  d^   =    Q(x)     , 
dx4 

y(0)  =  yd)  =  y-  (0)  =  y'  (1)  =  0  . 
From  (3.2.6)  and  (3.2.2)-(3.2.6),  it  follows  that 
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1 

(3.2.8)  e(x)  =    /  G^(x,z)  y(z)  dz 

0 

where 

z  (x-1)  ,  0  j<  z  <  X  _<  1 
G,(x,z)  = 

x(z-l)  ,  0  £  X  <  z  <  1 

is  the  Peano  kernel  for  linear  interpolation  used  in  the 
proof  of  Theorem  2.1. 

Similarly,  from  Birkhoff  and  Priver  (or  by  applica- 
tion of  the  Peano  theorem) ,  we  have 

1 

(3.2.9)  y(z)  =   /  G4  (z,t)  Q(t)  dt  , 

0 


where 


(3t2-2t^)z-^  +  3(t-2)t2z2 


+  3t2z  -  t^  ,   t  £  z 

^^4(2't)  =       2    3     3  2   2 

(3t^-2f^-l)z-^  +  3(t-l)^tz^  , 

t  >  z 

for  Oj<t_<l,   Oj<z_<l. 

Combining  (3.2.8)  and  (3.2.9)  ,  we  have 

1  1 

(3.2.10)   e(x)  =  /  G.  (x,z)      /  G4(z,t)  Q(t)  dtdz 
0  0 

=  /  G^(x,z)      /  G4(z,t)  u^^Nt)  dtdz 
0  0 

=   /    /  Gi(x,z)   G4(z,t)  u(^)(t)  dtdz 
0    0 

=   /    /  G^(x,z)   G4(z,t)  dz   u(^)(t)  dt 
0    0 
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-   /  G(x,t)   u(^) (t)  dt 


where 

1 

(3.2.11)  G(x,t)  =  /  G]_(x,z)  04(2, t)  dz  . 

0 

From  (3.2.11)  and  (2.2.5)  -  (2.2.8)  ,  it  follows  that 

(3.2.12)  g(2,0)  (j^^t)  =  G4(x,t) 
and 

(3.2.13)  g(P^2,0)  (j,^t)  =  G4(P'°)(x,t)  ,  p  =  0,  1,  2,  3  . 
Also,  as 

G4(z,t)    <.0,  0<z£l,      0<t<_l 

Gi(x,z)    10,  Oj<x<l,      0<z<l 

it   follows    that 

G(x,t)    >_  0  0<x<l,      0<t<l 

From  (3.2.10)  and  G(x,t)  >^  0,  we  have 

1  ,^ 

(3.2.14)  |e(x)|  <  /  G(x,t)  dt  maxg^^^^;,^  |  u  ^"^Mx)  |  . 

0  — 

In  fact, 

1  1 

(3.2.15)  /  G(x,t)dt  =      ///  G4(x,t)dt  dx  dx  +  ax  +  b 

0  0 

where  a  and  b  are  chosen  to  satisfy 

1  1 

(3.2.16)  /  G(0,t)dt  =  /  G(l,t)dt  =  0  . 

0  0 

We  know  from  Birkhoff  and  Priver  (or  Hermite)  that 
/  G.(x,t)dt  =  [-x2(l-x)2]  /  (4!)  . 


Then 


57 


/  G4(x,t)dt  dx  dx  =  _J. (  -  5x^  +  3x^  -  x^)  , 

0  6!       2 

and  to  satisfy  (3.2.16),  we  have  a  and  b  of  (3.2.15)  as 

a  =  1/2  ,      b  =  0  . 

6! 

Rearranging,  we  have 

(3.2.17)  /  G(x,t)dt  =  [  -  5x^^/2  +  3x^  -  x^  +  x/2  ]  /6! 
0 

=  [  x^(l-x)^  +  (l/2)x2(l-x)2  +  (l/2)x(l-x)  ]  /  6! 

Combining  (3.2.14)  and  (3.2.17)  ,  we  have  the  result  of  the 
theorem  for  p  =  0. 

From  3.2.10,  we  have 

(3.2.18)  |e(PMx)|</      |g(P)  (x,t)  I    dt   maxo<x<i|f  ^^^  (^)  I     • 

0  — 

From   3.2.11,    we   have 

X  1 

(3.2.19)  G<^'°Mx,t)    =/    yG4(y,t)dy+/     (y-1)  G4  (y ,  t)  dy    . 

0  X 

Therefore   as   G4(y,t)    <    0  0<y<l,      0<t<l 

X 


(3.2.20)        \G^^'^Ux,t)\    <I    y|G4(y,t)|dy 

0 


1 
+     /     (1-y)  iG4(y,t) |dy    . 

X 


As   before,    we   have 


/     |G4(y,t) |dt    =   y2(i-y)2    /    41     . 
0 
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Thus 

1  1  X 

(3.2.21)   /  |g(1'")  (x,t)  |dt  <  /  /  y  |G4(y,t)  |  . dy  dt 
0  0  0 

1  1 
+   /  /  (1-y)  |G4(y,t) I  dy  dt 

0  X 

X     1 
=   /  y  /  |G4(y,t) I  dt  dy 
0     0 

1  1 

+   /  (1-y)   /  iG4(y,t) I  dt  dy 
X        0 


=  /  [  y^(i-y)^  ]  /  4!  dy 
0 

+  J  [  (i-y)^y^  ]  /  4!  dy 

X 

=  [  1/60  -  (x3(l-x)^)/3  ]  /  4! 


=  fo,i(^) 

which  achieves  its  maximum  value  of  1/1440  for  x  =  0  or 
X  =  1.   We  note  also  that 

1/1440  =  1/(2  6!)  =  C]_ 

=  maxo<x<i|fo,o^^^  (^H  =  "^^^-0<x<l  1  ^  0  , 1  (^H  • 
Combining  (3.18)  and  (3.13),  we  have 

(3.2.22)   |e*P) (X) |  < 

1 
/     |G4(P-2'0)  (x-t)  |dt   maXo<x<i|u(^Nx)  1    , 
0  — 

for    p    =    2,    3,    4,    5. 

As  this  inequality  is  precisely  that  used  by  Birkhoff 

and  Priver  to  derive  the  functions  fg  2'    ^o  3'  -^0  4  ^'^^ 
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fQ  5,  the  theorem  follows  for  p  =  2,  3,  4,  5.   The  proof 
of  Theorem  3.2  is  very  similar  and  hence  omitted. 


CHAPTER  FOUR 
A  QUARTIC  SPLINE 


Introduction  and  Statement  of  Theorems 
Among  the  many  beautiful  properties  of  the  complete 
cubic  spline  is  the  fact  that  for  a  given  partition  and 
function  values,  the  cubic  spline  is  obtained  by  solving  a 
t r i d i agona 1 1 y  dominant  system  of  equations. 
Unfortunately,  when  one  uses  higher  order  complete  splines 
the  bandwidth  grows.  In  fact,  for  a  2m  times  continuous 
spline  of  order  2m+l,  the  bandwidth  of  the  system  of 
equations  is  2m+l.  Furthermore  the  diagonal  becomes  less 
dominant  as  k  increases. 

It  is  natural  then,  to  increase  the  order  of  the 
spline  but  preserve  bandwidth.  Ideally  we  would  hope  to 
increase  the  diagonal  dominance  and  order  of  convergence. 
In  this  chapter  we  introduce  a  quartic  C^  '  spline  which 
gives  0(h)  rate  approximation  to  a  c'^'  function.  The 
quartics  are  obtained  by  the  solution  of  a  tr idiagonal ly 
dominant  system.  As  desired,  it  is  more  diagonally 
dominant  than  the  system  associated  with  the  complete 
cubic  spline. 

The  main  result  of  this  chapter  will  be  to  give  an 
exact  error  bound  for  the  quartic  spline  discussed  here. 
We  first  give  the  definition. 
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Let  f  be  a  real-valued  function  defined  on  [a,b]. 

k 

Choose  a  partition  fXj^3j^_Q  such  that 

a  =  Xq  <  X]^  <  .  .  .  <  Xj^  =  b  . 
Let  Zj^  =  (Xj^_-]^  +  Xj^)/2,  be  the  midpoint  of  [Xj__-j_/  Xj_]  for 
i  =  1,  2,  .  .  .  ,  k  and  for  these  i  set  h^_-^   =  Xj_  -  Xj__-]_. 
Definition  4.1   Given  the  function  f  and  the  partition 
{Xj^}j^_Q,  we  define  a  quartic  spline  s  (x)  such  that 

(4.1.1)  s(x)  e  c2[a,b]  (^  P^Ui-i,    x^]  ,  i  =  1,  2,  .  ,  k  ; 
(where  P4tXj^_-j_,  x^]     denotes  the  functions  which  are 
quartics  when  restricted  to  ['^i-ir    Xj^]) 

(4.1.2)  s(x^)    =    f(x^)  for  i  =  0,  1,  .  .  ,  k 
s(Zj^)  =  f(zj^)  for  i  =  1,  2,  .  .  ,  k  ; 

and 

(4.1.3)  s'(a)  =  f'(a)  and  s'(b)  =  f'(b)  . 

Lemma  4.1   Let  f  be  a  real-valued  funtion  defined  on 
[a,  b]  and  let  fXj_}^_Q  be  a  partition  of  [a,b].   A 
quartic  spline  s  satisfies  Equations  (4.1.1)  and  (4.1.2) 
if  and  only  if  s  satisfies  the  tridiagonal  system  of 
equations  for  i  =  1,  2,  .  .  ,  k  -  1 
(4.1.4) 

-h^  s' (x^_^)    +  4(hi  +  h^_-^)    sMx^)  -  hi_]_  sMx^^^) 

=  -11   [(^i-l/h^)   -   (hi/hj__^)]   f(Xj^) 

+  16  [(^i-l/h^)  f(Zi+i)  -  (^i/h^.^)  f(Zi)] 
-  5  [(^i-i/h^)  f(x^^;L)  -  (^i/hi_i)  f(Xi_i)] 

where  h^^  =  ^^.  +  1   ~   ^i' 

We  will  give  the  proof  later. 
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Assuming  from  (4.1.2)  that  f(Xj_),  i  =  0,  1,  .  .  ,  k 

and  f(Zj^),  i  =  1,  2,  .  .  ,  k  are  known,  then  (4.1.4)  is  a 

system  of   k  -  1  equations  in  the  unknown  variables 

s'{Xj^),  i  =  1,  2,  .  .  ,  k.   If  we  impose  the  conditions  of 

(4.1.3)  that  sMxq)  =  f'(a)  and  s'(Xj^)  =  f'(b)  are  given, 

we  have  k-1  unknowns  and  the  k  -  1  diagonally  dominant 

equations  (4.1.4).   Lemma  4.1  thus  assures  us  that  s'(x^) 

can  be  uniquely  determined  for  given  conditions  (4.1.1)- 

(4.1.3).  As  will  be  shown  in  the  proof  of  the  lemma,  there 

is,  on  any  given  subinterval  [x^,  x^^2_]  ,  a  unique  quartic 

Sj^(x)  satisfying  the  five  conditions 

(4.1.5)      ^^(Xi)  =  f(x^)  ,    s^(Zi+^)  =  f(Zi+i)  f 

^i^^i  +  1^  =  f(Xi  +  l)  '   s'i^^i^  "  s'(xj_)  , 

s'i^^i+l)    =    s' (^i+l)     • 
Equations    (4.1.4)    are  derived  by   imposing  the  conditions 

that 

Sj^"(x^  +  )    =    ^i-i'(^i-)    for    i    =    1,    2,    .    .    ,    k   -    1    . 
For    i    =    1,    2,    .    .    ,    k,       s^(x)    is    thus    the    restriction   of 
the  spline  s  to    [x^,  x^^^-j^]. 

Summarizing,  unique  solution  of  (4.1.4)  implies  that 
s(x)  is  uniquely  defined  on  each  partition  subinterval 
[Xj^,  x^^-^],  i  =  0,  1,  .  .  ,  k-1,  which  is  to  say,  on  all 
of    [a,    b].      We   have   shown 

Corollary    4.1    The   quartic    spline    of    Definition    4.1    is,    for 
a   given   partition    Cx^}^_q    and    function    f,    unique. 

We  now  make  the  comparisons  with  the  complete  cubic 
spline     more     explicit.         The     system     of     equations 
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corresponding  to  (4.1.4)  for  the  complete  cubic  spline  has 
left-hand  side 

hj_s'(Xi_l)  +  2(hj_  +  h^.^)  s'(Xj_)  +  hi_;Ls'  (^i  +  l^  * 
In  comparison  (4.1.4)  is  twice  as  diagonally  dominant. 

To  interpolate  the  2k  +  1  function  values  f(x^)  and 
f(Zj^)  using  our  C^^^  quartic  required  solving  the 
tridiagonal  system  of  k  -  1  equations  (4.1.4).  As  the 
cubic  spline  must  match  derivative  and  second  derivative 
values  at  each  interior  function  value,  interpolation  of 
the  same  2k  +  1  function  values  by  the  C^^^  cubic  spline 
would  entail  solution  of  a  system  of  2k  -  1  equations.  In 
other  words,  the  matrix  equation  to  be  solved  for  the 
quartic  is  only  half  as  large  as  that  required  for  the 
cubic. 

We  can  now  state  the  main  theorem  of  this  chapter. 
Given  a  partition  Cxj^3^_q  of  [a,  b]  ,  denote 

h  =  maxo^i<k-ihi  =  niaxo<  i<k_i  (x^^i-x^)  . 
For  each  x  in  [a,b]  ,  there  exists  i  such  that 
0  <_  i  <_  k  -  1  and  x^^  <_  x  <_  Xj_  +  i-      We  set  t  =  (x  -  x^)/h^. 
Theorem  4.1   Let  f  e  C^^^[a,b]  and  let  Cxj^I^^q  ^^  ^ 
partition  of  [a,b].   Let  s  (x)  be  the  twice  continuously 
dif f erentiable  spline  corresponding  to  f  and  ^x^^Jj^^q, 
where  s  satisfies  (4.1.1)  -  (4.1.3)  .   Then 

(4.1.6)    |f(x)  -  s(x)i  <  |c(t)|  h^  max^<x<j3|f  (5)  (X)  I  /  5! 
where 

c(t)  =  [3t2{i-2t)  (l-t)2  +  t(l-2t)  (1-t) ]  /  6  . 
Define 


'0<t<l 
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c(t) 


=  (  /  1  -  _1  )  (1/5  +  2/30  )  /  6  . 

4   /To 
It  follows  that 

(4.1.7)  |f  (X)  -  s(x)  I  <  Cq  h^  max^<j^<i3|f  ^^)  (x)  |  /  5!  . 
Furthermore,  neither  |c(t)]  nor  c^  can  be  improved,  as  we 
can  show  by  letting  f  =  x^/5!  and  letting  k  become 
arbitrarily  large  for  an  equally  spaced  partition.  An 
approximate  decimal  expression  for  c^  is  .0244482  and 
Cq/5!  is  approximately  .000203818  . 

We  will  also  show 

(4.1.8)  |f  '  (x^)  -  s'  (Xj_)  I 

<  h"^  max^<x<bl^^^^  (^)  I  /  ^• 
and  that  this  estimate  is  exact. 

Related  to  Theorem  4.1  is  the  following  conjecture. 
Conjecture  4.1   Let  f  6  C^^^[a,b]  and  let  Cxj_}^^q  be  a 
partition  of  [a,b].   Let  s(x)  be  the  twice  continuously 
dif f erentiable  spline  corresponding  to  f  and  Cxj^}|^_q, 
where  s  satisfies  (4.1.1)  -  (4.1.3  )  .   Then 

(4.1.9)  |f'(x)  -  s'(x)|  <  h^max^^j^^j^lf  (^)  (x)  I  /  6!  . 

If  Conjecture  4.1  holds,  then  the  constant  1/6!  can  not  be 
improved.  This  conjecture  has  been  verified  numerically. 
Remark  4.1  Given  f  6  C  [a,b]  and  a  partition  fXj^3^_Q  of 
[a,b],  let  s  be  the  quartic  C^^^  spline  satisfying 
(4.1.1)  -  (4.1.3)  .  Then  the  supremum  norm  Mf^"""^  -  s^-*-^!! 
is    of    order    h^~^    max^^-^^^j^  |  f  ^^^  (x)  |  ,        i    =    0,    1,    2. 
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Theorem  4.1  demonstrates  that  the  quartic  C^  '  spline 
gives  the  best  possible  order  of  approximation  to 
functions  from  the  smooth  class  C^^' .  We  next  discuss 
interpolation    to    the    much    less    smooth    class    of    functions 

which  are   merely  continuous   on    [a,b].      As   f '  (a)    and   f'(b) 

(2) 
are    not   necessarily   defined,     we    consider    the    quartic    C^    ' 

spline     satisfying     (4.1.1)     and     (4.1.2)     with    boundary 

conditions 

(4.1.10)  s' (a)    =    s' (b)    =    0    . 

Denote  (^(f/h)  =:  sup  i  jj_y  i  <j^  |  f  (x)  -  f(y)l. 
Theorem  4.2   Let  f  6  C[a,b].   If  Cxj_}^^o  ^^  the  partition 
of  equally  spaced  knots,  then  for  Xj^  <_  x  <_  Zj_  +  i    =    (^^    + 
Xj_  +  ]_)/2  and  t  =  (x  -  x^) /h^,     i  =  0 ,  1 ,  .  .  ,  k  -  1 ,  we 
have 

(4.1.11)  |f  (x)  -  s(x)  I  <    c(t)  aj(f,h)  ;  0  <  t  <  1/2 
and  for  Zj^^-j^  i.  x  <_  x^^-j^,  or  1/2  <.  t  _<  1 

(4.1.12)  |f (x)  -  s(x) I  <  c(l-t)  w(f,h) 

where  c(t)  =  [1  +  (13/3)t  -  Bt^  -  (58/3)t^  +  let"^}  . 
Note  that  maxg^^^-j^ /2C  (t)  is  approximately  equal  to 
1.6572. 

The  bound  of  the  preceding  theorem  is  only  valid  for 
equally  spaced  knots.  For  arbitrary  partitions  we  can  not 
give  a  bound  of  this  same  form.   However,  if 

m  =:  maxo<i<_k-ihi  /  niino<_i<}c-ihi  ' 
we  have  the  following  theorem. 
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Theorem  4.1.3   Let  f(x)  e  C[a,b],  and  let  s  be  the  C 


(2: 


quartic  spline  satisfying  (4.1.1),  (4.1.2),  and  (4.1.9). 
Then  for  X-;  <.  x  <_  z^  +  i/  and  i  =  0 ,  1 ,  2 ,  .  .  ,  k  -  1 , 
(i.e.,  for  0  <  t  <_  1/2  with  t  =  (x  -  x^) /h^) 

(4.1.13)  |f(x)  -  s(x)|  <  C;l<^)  oa(f,h)  , 
and  for  z^^-j^  £  ^  £  ^i  +  i'  i«e.  1/2  <.  t  _<  1, 

(4.1.14)  |f(x)  -  s(x)|  <  C3_(l-t)  u)(f,h) 
where 

c^(t)  =  [1  +  lOt^  -  28t^  +  let"*] 

+  (8/3) [m^  +  m]  [t (l-2t) (1-t) ]  . 
Theorems  4.2  and  4.3  indicate  that  for  suitable 
partitions  the  quartic  C   spline  can  provide  acceptable 
approximations  to  functions  which  are  merely  continuous  on 
[a,b]  . 

Proof  of  Lemma  4.1 
We  first  give  an  expression  for  the  unique  quartic 
matching  function  and  derivative  values  at  endpoints  and 
function  values  at  the  midpoint.  Specifically,  let  f  be  a 
real-valued  function  defined  on  [0,1],  and  dif f erentiable 
at  0  and  1.  Let 
(4.2.1) 

P^(x)    =  1  -  llx^  +  18x^  -  8x'*  =  (l-2x)  (l-x)2(l  +  4x)  , 
P2(x)  =  IGx^  -  32x^  +  16x^  =  IGx^d-x)^  , 
P3(x)  =  -Sx^  +  14x^  -  8x^  =  -(l-2x)x^[l+4(l-x) ]  , 
P4(x)  =  X  -  4x2  ^  5j^3  _  2x'^  =  x  (l-2x)  (1-x)  ^  , 
P5(x)  =  x^  -  3x^  +  2x'^  =  x2(l-2x)  (1-x)  . 
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Then 

(4.2.2)  L[f,x]  =  P^(x)f(0)  +  P2(x)f(l/2)  +  P3(x)f(l) 

+  P4(x)f ' (0)  +  P5(x)f ■ (1) 
is  the  unique  quartic  satisfying 

(4.2.3)  L[f,0]  =  f(0)  ,  L[f,l/2]  =  f(l/2)  , 

L[f,l]  =  f(l)  ,  L'[f,0]  =  f'(0)  ,  L'[f,l]  =  f  (1)  . 
L  is  a  linear  functional  and  a  projection.  If  f  is  a 
polynomial  of  degree  four  or  less,  then  L[f,x]  =  f(x).  In 
the  future  calculations  we  will  need  the  following  facts 
about  the  quartics  P-j_, 


(4.2.4) 


(0)  =  -22 

(0)  =  32 

(0)  =  -10 

(0)  =  -8 

(0)  =  2 


(1)  =  -10 
(1)  =  32 
(1)  =  -22 
(1)  =  -2 
(10)  =  8 


Let  z^^-]_  =  (x^  +  Xj_  +  -l)/2.  On  the  interval  [x^,  ^i  +  1^' 
the  unique  quartic  Lj_[f,x]  interpolating  f(x^),  f  (x^)  , 
f(z^_^,-^),  f(x^^-^),  and  f  {x^,^.-^)  can  be  expressed  in  terms 
of  Pj_.  In  fact,  let  t  =  (x  -  x^) /h^  where  hj_  =  ^i  +  i  -  ^j_- 
Then 
(4.2.5)    Lj_[f,x]  =  f(x^)    P^(t)  +  f(Zi  +  i)  ^2^^^ 

+  f(Xi  +  ^)  P3(t)  +  h^  f'(Xj_)  P4(t) 

+  h^  f  •  (x^^3_)  P5(t)  . 
Let  s  be  the  quartic  spline  of  Definition  4.1  corres- 
ponding to  f  and  the  given  partition.   Then  the  restric- 
tion s-  (x)  of  s  to  [x  j_  ,Xj_^-]_  ]  is  a  quartic.   Hence 
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L^[s,x]  =  s^(x).  Using  the  facts  that  s(Xj^)  =  f(Xj^), 
s(Zj_^-l)  =  f(zj_^j_),  and  flXj^^.-^)  =  s(xj_^^)  in  (4.2.5)  we 
have 

(4.2.6)  s^(x)    =  f(x^)  P^(t)  +  f(Zi+i)  P2<t) 

+  f  (x^^3^)  P3(t)  +  h^  s'  (x^)  P4(t) 
+  h^  s'(x^^^)  P5(t)  ,  t  =  (x-Xi)/h^  . 
In  order  that  s  be  twice  continuously  dif f erentiable, 
we  must  satisfy 

(4.2.7)  s^' '  (x^+)  =  s^_i' '  (x^-) 

where  s^  is  the  restriction  of  s  to  [x^,  ^i  +  i]  and  s^_-^  is 
the  restriction  of  s  to  [Xj^_-|_  ,Xj_  ]  .  Differentiating 
(4.2.6)  twice  we  have 

(4.2.8)  s"(x^  +  )  =  -J^  C  f(Xi)P^"(0)  +  f  (z^^^)P2"  (0) 

+  f(Xj^^^)  P3"(0)  +  h^  s'(Xj_)  P4"(0) 

+  h^  s-  (x^^;^)  ^5*  '  (0)  5  • 
Similarly,  from  rewriting  (4.2.6)  for  the  interval 

['^2.-1'    ^i-'  '  ^®  have 

(4.2.9)  s"(x^-)  =  — L^  [  f  (Xj^_^)P^' •  (1)  +  f(z^)P2"(l) 

^i-l 

+    f (x^)    P3' • (1) 
+   hi_i    s-  (x^.^)    P4"  (1) 
+   h^_^   s' (x^)    P5' •(!)}. 
Setting     s"(Xj_  +  )      =     s"(Xj^-)     by     eguating     (4.2.8)     and 
(4.2.9)     and     using     Pj^"(0)     and     Pj_"(l)     from     (4.2.4),      we 
have 
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(4.2.10)  [-22    f(Xj_)     +    32    f(Zi  +  i)     "    10    fiXj^  +  j^) 

-    8    h^    s'(x^)     +    2    h^    s'(x^^^)     }    /    h^^ 
=    C-10    f(x^_;L)     ■"    32    f(z^)    -    22    f(Xj_) 

-  2    hj__]_    s'(x^_^)     +    8    hj__3^    s'(x^)     ]    /    h^_-^^    . 
Factoring   two,    multiplying   by   hj^h^_-|_,    and   putting   the 

known  function  values  on  the  right  hand  side,  we  have 

(4.2.11)  -h^s'(x^_^)    +    4(h^+h^_^)    s'(Xi)    -   ^^-i^'  i^i+i) 

=   -11    [(^i-i/h^)    -    (^i/h^_^)]    f(Xj_) 
+    16    [(^^i-l/h^)    f(Zi  +  i)    -    (^i/h^_^)    f(z^)] 

-  5  [(^i-l/h^)  f(x^  +  ;L)  "  (^i/hi_i)  f(Xi_i)]  , 
which  is  the  desired  system  of  equations  (4.1.4).    Having 
established  Lemma  4.1,  we  next  turn  to  a  proof  of  Theorem 
4.1. 

Proof  of  Theorem  4.1 

Our  method  of  proof  is  to  establish  a  pointwise 
bound.  As  in  the  proof  of  Lemma  4.1,  let  Lj^[f,x]  be  the 
unique  quartic  agreeing  with  f(xj_),  f(x^^.]_),  f(Zj_^-[_), 
f'(Xj_),  and  f'(Xj_^-[_),  and  let  s  be  the  twice  continuous 
quartic  spline  corresponding  to  f  and  Equations  (4.1.1)  to 
(4.1.3)  on  the  partition  fXj_}^^Q.  Then  for  Xj_  <_  x  <_  x^  +  ir 
we  have 
(4.3.1)    lf(x)  -  s(x)|  <  |f(x)  -  L^[f,x]| 

+  |L^[f,x]  -  s(x) I  . 

Assume  that  f  6    C^^^  [a,b].      By  a  proof  attributed  to 
Cauchy,  we   know   that 
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(4.3.2)  lf(x)  -  L^[f,x]|  <  (h^5/5!)  | t^ (1/2-t) (1-t) ^ |  U 
where  t  =  (x-x^j/h^^  and  U  is  the  maximum  of  |f^^Mx)|  on 
[x-,x.^^].   Equation  (4.1.9)  gives  a  pointwise  bound  for 
|f  (X)  -  Lj_[f,x]  |. 

Let  i  be  arbitrary  and  x^  <_  x  <_  x^^^^.   We  next  turn 
our  attention  to  deriving  a  similar  bound  for 
|L^[f,x]  -  s(x)|  =  |L^[f,x]  -  Sj^(x)|.   Subtracting  (4.2.6) 
from  (4.2.5)  gives 

(4.3.3)  Lj_[f,x]  -  s^(x)  =  h^  [f'(xj^)  -  s'(Xj_)]  P4(t) 

+  h^  [fMx^^l)  -  sMx^^^)]  P5(t)  . 
Denoting 

(4.3.4)  e'  (x^)    =   f '  ix^)    -   s'  (x^^)  , 
then  we  have  from  (4.2.5), 

(4.3.5)  |L^[f  ,x]  -  Sj_(x)  I  < 

h^  max(  ieMx^)  I  ,  |e'(x^^^)|3  [  |  P^  (t)  |  +  |  P5  (t)  |  }  . 
As  P4(t)  =  t(l-2t)  (1-t)^  and  P5(t)  =  t^  (l-2t)  (1-t) 
are  both  positive  for  0  <_  t  <_  1/2  and  both  negative  for 
1/2  <  t  <  1,  lP4(t)  I  +  |P5(t)  I  =  |P4(t)  +  P5(t)  I  for 
0  <_  t  <_  1.   Then  for  Xj_  <_  x  <_  x^^-|_,  we  have 

(4.3.6)  |Lj_[f  ,x]-s(x)  I  < 

h^  maxC|e'  (x^)  I  ,|e'  (Xj_^l)  |}|t(l-2t)  (1-t)  I  . 
Redefine  L  so  that  its  restriction  to  [Xj_,  Xj_^.j_]  is  L^^  for 
each  i,  i  =  0,  1,  .  .  ,k-l.  Choose  i  so  that  |e'(Xj_)|  is 
maximal.  We  then  have  for  all  a  <  x  <  b. 
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(4.3.7)  |L[f,x]  -  s(x)|  <  h  |e'(x^)|  1 1  ( 1/2-t)  (1-t)  | 
where  h  =  maxQ^^^]^_-i_h^  is  the  maximal  subinterval  length 
and  where  on  each  subinterval  [Xj,  Xj^^.^'  "^  <.  J  1  ^  ~  ^' 
we  define   t  =  (x  -  Xj)/hj. 

The  next  task  is  to  bound  |e'(xj_)  |.   From  both  sides 
of  (4.1.4)  we  subtract 

-h^  f'(Xi_l)  +  4(h^   +   hi_j_)f'(x^)  -  h^_^f'  ix^^^)     , 
thereby  defining  a  functional  BQ(f) 

(4.3.8)  hj^e'(Xi_i)  -  4  (h^+h^.^)  eMx^)  +  h^.^e' (Xi^]_) 
=  h^  f'(Xi_i)  -  4(hi  +  h^_i)f  (x^)  +  hi_if'(x.^i) 

-11  [(^i-l/h^)  -  (^i/h^.^)]  f(Xi) 
+  16  [(^i-l/h^)  f(Zi+i)  -  (^i/hi_i)  f(Zi)] 
-  5  [(^i-l/h^)  f(x^^i)  -  (^i/h^_]_)  f(x^_i)] 
=:  Bo(f)  . 
The  linear  functional  BQ(f)  is  identically  equal  to 
zero  when  f  is  a  polynomial  of  degree  four  or  less,  as  can 
be  directly  verified.   (The  arithmetic  of  verification  is 
simplest  if  one  takes  x^.-^  =  -'^]_-\r   Xj_  =  0,  x^^-|_  =  hj_  and 
checks  the  monomials  1,  x,  x^,  x^,  and  x  ). 

We  have  chosen  i  so  that  |e'(Xj_)|  attains  its  maximum 
value.   As 

4(h^  +  hj^_3^)e'  (x^) 

=  -BQ(f)  +  h^eMx^.^)  +  hj^_;L^'<^i  +  l^' 
it  follows  that 

|4(hi  +  h^_^) I  <  |Bo(f)|  +  |h^e'(x^_^)| 

+  |h^_^e' (x^) 

<  |Bo(f)  I  +  1  (hi  +  h^.^L^^'  (^i^  I  • 
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Hence 

and 

(4.3.9)  le'(Xi)|  <  |Bo(f)|  /  [3  (hj^+h^.^)  ]  . 

As  BQ(f)  is  a  linear  functional  which  is  zero  for 
polynomials  of  degree  four  or  less,  we  can  apply  the  Peano 
theorem  to  get 

(4.3.10)  BQ(f)  =   /  ^^■^Bo[(x-y)^^]  f^^Ny)  dy  /  4 !  . 

^i-1 
From  (4.3.10)   follows 

(4.3.11)  |BQ(f)|  <    /  ^^•^|BQ[(x-y)^4]  I  dy  U^/  4!  . 

^i-1 
where  Uj^  is  the  maximum  of  |f^  '|  on  [Xj^_-|  x^^-j^]. 

For  x^_-]_  <.  y  <.  ^i  +  i'  BQ[(x-y)^^]  takes  the  form 

(4.3.12)  Bo[(x-y)^4]  = 

-  16(h.+h._^)  (Xi-y)^3  +  4h^_i  (x^^l-y)3 

-  11  [(^i-l/h^)  -  (hi/h^_^)]  (x^  -  y)+^ 

+    16    [(^i-i/h^)  (z.^^-y)^4    -    (l^i/hi_i)  (z^-y)^^ 
-   5    (^i-l/h^)     (x^^^   -  y)^    . 
In    order    to    evaluate    the    integral    of    (4.3.11)    we    need 
to    know    the    sign    behavior    of    Bq  [  (x-y)  ^'*  ]  .       We    rewrite 

(4.3.12)  in    a    form    which    shows    its    symmetry   about   x^. 

(4.3.13)  Bo[(x-y)^^]    = 

(^i/h^_^)     [-5(x^-y)    +   h._^]     [(x^   -  y)    -  hi_i]3    , 


for   x^ 


-1    <  y  1  ^i 
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for  Zj^  £  y  _<  Xj^ 
(^i-l/h^)  {Xj_-y)2  [ll{Xj_-y)2  +  16h^{x^-y) 

for  x^  <  y  <  Zj_^^ 
(^i-l/h^)  [-5(x^-y)  -  h^]  [(Xj^-y)  +  h^]^  , 

for  z^^^  <  y  <  x^^i  . 
As  the  expression  (4.3.13)  has  factors  which  are  at 
most  quadratic  it  is  fairly  easy  to  to  determine  to 
determine  the  sign  of  Bq  [  (x-z)  ^.'^  ]  .  In  fact,  Bq[(x-z)+  ] 
is  nonnegative  for  ^_-^  1  y  1  ^i  +  l-  Evaluation  of  (4.3.11) 
is  then  straightforward.   The  term  by  term  integration  of 

(4.3.13)  gives 

(4.3.14)  /''^^^|Bo[(x-y)/]  I  dy  =  h^h^.^^  ^^i-l^  ^  h^^J/lO  . 
^i-1 

From  Equation  (4.3.11)  we  conclude  that 

(4.3.15)  lBQ(f)|  <   U^  h^hi_i[hi_i^  +  h^^]  /  [2(51)]  . 
From  (4.3.9)  it  is  then  evident  that 

(4.3.16)  |e'  (Xj)  I  < 

"i  ^i'^i-l^'^i-l^  ^  *^i^^  /[(6!)  (h^+h^_^)] 
for  j  =  1,  2,  .  .  ,  k-1.   As 

[h^_;L^  +  '^i^^  /  [hj^+h^_]^]  <  max  ih^^ ,    h^^.^^}  , 
and  as 

Uj_  <  U  , 
it  follows  that 
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(4.3.17)  max    |e'(Xj)|    <   maxChj^^,    h^_^^'i    U/(6!)     . 
This     is     the     desired    bound     on     |e'(Xj^)  |  . 

Applying   it  in    (4.3.7)    we  have 

(4.3.18)  |L[f,x]    -    s(x)|    <    h^    |t(l-2t) (1-t) I    U/     (6!)     . 
From    (4.3.2)    follows 

(4.3.19)  |f (X)    -    L[f,x] I    <    h^     |t2(t-l/2) (l-t)2|    U/    5! 
where   L   restricted   to    [x^,    Xj^^-]^]    is   defined   as   Lj^[f,x]    and 
where  h  is  the  maximum  of  hj  . 

We   can   now   combine   the   bounds   on    |f(x)    -   L[f,x]  |    and 
|L[f,x]    -    s(x)|.      Prom    (4.3.19)    and    (4.3.18),   we   have 

(4.3.20)  |f(x)    -    s(x)|    <    h^     |c(t)|    U    /    5! 
where 

lc(t)|    =    |3t2(i-2t) ((l-t)2|    +    |t(l-2t) (1-t) I    /    6 
=    |3t2(l-2t) (1-t)^    +    t(l-2t) (1-t) I     /    6 
and 

c(t)       =    [3t(l-t)    +    1]     [t(l-2t) (1-t)]    /    6    • 
Then 

(4.3.21)  c^   =   maXQ<^<^|c(t) |    • 

To   verify    (4.3.21),    note    that 

(4.3.22)  6c'(t)     =    -30t2(t-l)2    +    1 

=    -30     [(t   -    1/2)     +    l/2]2     [(t    -    1/2)    -    l/2]2    +    1 
=    -30     [ (t   -    1/2)2    _    i/4]2    ^    1    . 
For    0    <    t   <    1,    the   roots   of   c" (t)    are 


(4.3.23)  t    =    1/2    +     /1/4    -    1//30    . 
Evaluating   c(t)    at   the   roots   of   c'(t),    we   get 

(4.3.24)  Cq    =     (  /1/4    -    1//I0    )      (1/5    +    2//30'    )     . 
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We  have  shown  the  so-called  "direct"  part  of  the 
proof,  that  Equation  (4.1.7)  holds  for  c^.  It  remains  to 
be  shown  that  the  theorem  holds  for  no  smaller  c^. 

In  fact,  given  c  <  c^,  we  can  produce  a  function  f 
and  a  partition  Cxj^}^^q  of  [-1,  1]  such  that 

(4.3.25)  max_2^<j^<]^  I  f(x)  -  s(x)|  > 

c  h^  max_i<jj<i|f(^^  (X)  |/5!  . 

Often,  when  polynomial  interpolation  of  degree  n  is 
considered,  the  worst  error  is  attained  by  a  polynomial  of 
degree  n  +  1.  As  s  is  a  quartic  spline,  it  is  natural  to 
try  f (x)  =  x^/5!  as  a  possible  worst  function.  A  particu- 
larly pleasant  feature  of  the  trial  worst  function  f  is 
that  it  has  fifth  derivative  identically  equal  to  one. 

For  X-  _<  X  _<  ^i  +  ir    we  have  by  the  Cauchy  formula 

(4.3.26)  x^/5!  -  L^[x^/5!,x] 

=  h^5[t2(t-l/2) (t-l)2]  /  5!  . 
Furthermore,  for  equally  spaced  knots  Xj__.j^,  x^^,  ^i  +  ir    we 
can  calculate 

(4.3.27)  Bq(x5/5!)  =  h^^ /    5!  . 

If  e'(x^_^)  =  e'(x^)    =  e'(x^^-^),    we  have  from  (4.3.8) 

(4.3.28)  e'(x^)  =  -Bq(x^)  /  6  =  -h^^/    6!  . 
Equation  (4.3.3)  then  becomes  for  f(x)  =  x-'/5! 

(4.3.29)  h^[f,x]    -    s(x)  =  -h^   h^^    i:P4(t)  +  P5(t))/6! 

=  h^5  [t(2t-l)  (l-t)]/6!  . 
Combining  Equations  (4.3.26)  and  (4.3.29),  we  have, 
for  x>  <_  x  _<  ^^  +  1' 
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(4.3.30)  f(x)  -  s(x)  =  h^^    Ct(t-l/2) (l-t)/3 

+  t2(t-l/2) ( (l-t)^}  /  51 
As  (4.3.30)  gives,  after  taking  its  absolute  value, 
precisely  our  pointwise  bound  |c(t)  |  of  (4.3.20),  we  will 
have  attained  c^,  provided  only  that  h^   =   h,  and  as  men- 
tioned above, 

(4.3.31)  e'(x^)  =  e'(x^^^)  =  e'(Xj__^)  =  -h^/6! 

In  order  that  h^^  =  h,  we  take  the  knots  to  be  equally 
spaced.  Attaining  (4.3.31)  is  not  so  easy.  In  fact  it  is 
attained  only  in  the  limit.  The  difficulty  is  the  boundary 
conditions  e'(xQ)  =  e'(X]^)  =  0.  We  can  show,  however, 
that  as  one  moves  many  subintervals  away  from  the 
boundaries,  e'(Xj^)  goes  to  -h'^/6!. 

Explicitly,  let  fxj^j^^g  be  the  partition  dividing 
[-1,1]  into  k  equal  subintervals;  in  this  case,  h  =  h-  = 
2/k  .   For  i  =  1,  2,  .  .  .,  k  -  1,  and  f  =  x^/5!,  we  have 
BQ(f)  defined  on  [Xj__-,_,  x^^^-j^]  and 

(4.3.32)  Bo(f)/h  =  h'^/B!  =  e'  {x^_^)    -    8e'(Xj^)  +  e'  (x^^-^) 
We  wish  to  apply  (4.3.32)  inductively  to  move  away 

from  the  end  conditions  e'(-l)  =  e'(l)  =  0.  In  order  to 
do  so  we  must  establish  that  e'(x^)  <_  0  for  0  _<  i  <_  k.  We 
reason  by  contradiction. 

Let  1  1  i  <.  k  -  1.   Suppose  e'(x^)  >  0.   Then 
e'  (x^_^)    +   e'  ix^  +  i)       > 

e' {x^_^)    -    8  e' (x^)  +  e' (x^^^) 
>   h'^/B! 
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Hence 

,4 


max[  |e'(x^_i)|,  |e'(x^^i)|  3>h^/[2(5!)]  , 
contradicting  the  fact  (4.3.17)  that 

h'^/ei   >  max  [  |e'(Xj^_^)|,  |e'(x^^^)|  ]  . 
We  have  shown  by  assuming  the  contrary  that 
e'(x^)     £  0  for   i  =  1,  2,  .  .  ,  k  -  1  . 
Condition  (4.1.3)  is  that  eMxg)  and  e'(xj^)  are  zero.   Thus 
(4.3.33)   e'(Xj^)  _<  0  for  i  =  0,  1  ,  .  .  ,  k  . 

Applying  (4.3.32)  again  we  have  for  i=l,  2,  ..,k-l 
8e'(Xj_)  =  -h^/S!  +  e'(Xi_^)  +  eMx^^^)  . 


As 


e'  (x^_-^)  ,    e'(Xj^^-[^)  _<  0,  this  implies  that 


ie' (x^)  <  -h^/5! 


and 


(4.3.34)  e' (x^)  <  -h^/ [8(5!) ]  . 

Similarly,  for  i=2,  3,  ..,k-2,  we  have 
8e'(x^)  =  -h^/5!  +  e'(Xi_i)  +  eMx^^^)  , 
and  hence  by  (4.3.34), 

eMx^)  <  -h^/S!  -  h^/[8(5!)]  -  h^/[8(5!)] 
=   -(1  +  1/4)  h^  /[8(5!)]  . 
Inductively,  for  i=jtoi=k-j,we  will  have 
e'(x^)  <  -{1  +  1/4  +  1/4^  +  .  .  .  +  l/4^"l3h^/ [8(5!) ]  . 
The  harmonic  series  (1  +  1/4  +  1/4^  +  .  .  .  )  is 
equal  to  1/(1  -  1/4)  or  4/3  .   Thus,  in  the  limit  as  i,  k, 
and  j  go  to  infinity,  we  have 

(4.3.35)  e'  (x^)  <  -(1/8)  (4/3)h'^/5!  =  -h^/6!  . 

We    already    know    from     (4.3.17)     that     |e'(Xj^)|     <    h^/6!.       Thus 
for   k   >    2j    +    1,    and   k-j>i>jasj    goes    to    infinity. 
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we  have 

(4.3.36)   e'(Xj^)  goes  to  -h^ / 6 1     . 

In  the  sense  of  (4.3.36),  (4.3.31)  is  satisfied. 
Then,  as  e'(xj_)  goes  to  -h'*/6!,  x^/5!  -  s(x)  goes 
uniformly  to  the  expression  of  (4.3.30)  and  (4.3.20)  with 
h  =  hj^.  It  follows  that  the  expression  of  (4.3.20) 
cannot  be  improved  further.   In  fact  we  have  shown  that 

|c(t)|  offers  a  pointwise  exact  bound,  and  its  maximum  c^ 
is  the  exact  norm  bound. 

Proof  of  Theorem  4.2 
We  know  from  (4.3.1)  that  for  Xj^  <_  x  <_  Xj_^-|^, 

(4.4.1)  s(x)  -  f(x)  =  P-L(t)  f(Xi)  +  P2(t)  f(Zi  +  ;L^ 

+  P3(t)  f(Xi^i)  +  h^  P4(t)  s'(x^) 
+  h^  P5(t)  s- (x^^^)  -  f (X)  . 
It  is  easily  verified  that  P-^it)    +   P2(t)  +  P3(t)  =  1. 
Thus 

(4.4.2)  s(x)  -  f(x)  =  P-^it)     [f(x^)  -  f(x)] 

+  P2(t)  [f  (z^^;^)  -  f  (X)] 
+  P3(t)  [f(x^^^)  -  f(x)] 
+  h^  s' (x^)  P4(t) 
+  h^  s'  ix^^-^)    P^(t)     . 
Each  of  the  first  three  terms  on  the  right  hand  side 
can  be  bounded  in  absolute  value  by    (f,h).   We  also  must 
bound  the  last  two  terms.  For  equally  spaced  knots  h  =  hj_ 
=  hjj^_-|_;   equation  (4.1.4)  reduces  to 
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(4.4.3)  -h  s'(x^_^)  +  8h  s'  (x^)    -   h   s'(x^^-|_) 

=  16  [f(Zi^i  -  flz^)] 

-  5  [f(Xi^3^)  -  f{x^_^)]     . 
Assume  s' (Xj^)  is  maximal  in  absolute  value.   Then 
|6h  s'  (x^)  I  <  16  |f  (Zi  +  i)  -  f (z^)]  I 

+  5  |f(Xi^i)  -  f(x.) I 
+  5  |f  (x^)  -  f  (Xj__i)  I 
<  26  w (f ,h)  , 
and  hence 

(4.4.4)  |h  s'(x^)|  <  (13/3)  aj(f,h)  . 
Combining  (4.4.2)  and  (4.4.4) ,  we  have 

(4.4.5)  |s(x)-f(x)|  <  [|Pi(t)|  +  |P2(t)|  +  |P3(t)| 

+  (13/3)  |P4(t)  I  +  (13/3)  lP4(t)  I  3  aj(f,h)  . 
For  0  <  t  <  1/2, 

(4.4.6)  |Pi(t) I  +  |P2(t) I  +  |P3(t) I 

+  (13/3)  !P4(t)  !  +  (13/3)  |P4(t)  I 
=  P^(t)  +  P2(t)  -  P3(t)  + 

(13/3)P4(t)  +  (13/3)P5(t) 
=  1  +  (13/3)t  -  3t2  -  (58/3)t^  +  16  t^  . 
We  have  shown  the  theorem  for  0  <_  t  <_  1/2.  The 
argument  for  1/2  <.  t  _<  1  is  symmetric. 

Proof  of  Theorem  4.3  We  are  considering  now  the  case  in 
which  knots  are  no  longer  assumed  to  be  equal.  We  assume 
that  the  ratio  of  the  longest  subinterval  to  the  shortest 
is  less  than  m.  Equation  (4.4.2)  still  applies.  Again  we 
choose  i  so  that  |s'(x^)l  is  maximal.  From  (4.1.4)  we  now 
have 
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(4.4.7)     -h^sMx^.^)  +  4(h^+h^_^)    s'(Xj_)  -  i^i-is'  (Xj_  +  ;L> 
=  -11  [('^i-l/h^)  -  (^i/hj__i)]  f(x^) 
+  16  [(^i-l/h^)  f(Zi+i)  -  (^i/hi_i)  flz^)] 
-  5  [(hi-l/hj_)  f(x^^^)  -  (hi/h^_^)  f(Xj__;L)l 
=  ll[^i/hi_i]  [f(Xi)  -  f(Zi)] 

+  5[^i/hi_i]  [-£(2^)  +  f(Xi_i)] 
+  ll[^i-l/hj^]  [f(Zi^i)  -  ffx^)] 


+  SL^^i-l/h^]  [f(Zi^i)  -  f{Xi^i)] 


Then 


3(h^  +  h^.^)  |s'  (x^)  I  <  16  ^i/hi_i  aj(f,  h^_^/2) 

+  16  ^i-l/hj_  oj  (f  ,h^/2) 
and 
ls'(x^)|  <  (16/3)  [^i/h^.^L  ■"  ^i-l/hi]   oj(f,h) 


hi  +  hi_i 

where  h  =  ni^Xg^j^^j,_-Lhj^ . 

Then  for  any  given  j,    0  <_  j  £  k-1  ,  and 

m  =  maxChj_}/ Cminj_} ,  i  =  0,  1,  .  .  .  ,  k  -  1,  we  have 

(4.4.8)    maxC|hj  s'(Xj)|,  |hj  s'(Xj^.3_)|] 

<^  (16/3)     h       (m  +  1/m)  oj(f,h) 
2min[hj^} 

<  (8/3)  (m^  +  m)  oj(f,h)  . 

Substituting  (4.4.8)  into  (4.4.2)  yields  the  result 

of  the  theorem. 


CHAPTER  FIVE 
IMPROVED  ERROR  BOUNDS  FOR  THE  PARABOLIC  SPLINE 


Introduction  and  Statement  of  Theorems 
The  quartic  splines  of  Chapter  Four  share  and  improve 
many  of  the  properties  of  the  complete  cubic  spline.  To 
insure  a  good  approximation  to  a  given  continuous 
function,  we  must  make  the  largest  subinterval  of  a 
partition  small.  Unfortunately,  we  must  also  pose  some 
additional  restrictions  on  the  partition.  For  instance, 
in  Theorem  4.3,  the  norm  of  the  error  depends  not  only  on 
the  length  h  of  the  largest  subinterval  but  also  on  the 
ratio  m  of  the  largest  to  smallest  length  subinterval. 
Similar  additional  restrictions  must  be  made  for  the  cubic 
spline. 

In  this  chapter,  we  will  discuss  a  spline  operator 
for  which  the  norm  of  the  approximation  error  goes  to  zero 
with  the  length  of  the  largest  subinterval,  for  any  par- 
tition and  any  continuous  periodic  function.  This  spline 
is  the  piecewise  parabolic  spline  introduced  by  Marsden 
and  discussed  in  Chapter  One.  Its  properties  are 
summarized  in  Equations  (1.6.1)  to  (1.6.7). 

As  Marsden  points  out,  many  of  the  bounds  he  gives 
can  be  sharpened.  The  main  result  of  this  chapter  will  be 
to  accomplish  this  sharpening.   While  many  of  the  bounds 
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given  here  may  still  not  be  exact,  at  least  one  of  them 
is,  and  in  fact  is  even  pointwise  exact.  In  other  cases  we 
can  reduce  the  known  bounds  by  a  factor  of  more  than  two. 

The  results  given  here  thus  enable  one  to  compare  the 
error  of  the  Marsden  spline  to  the  error  of  other  spline 
interpolation  processes.  Specifically,  future  work  on  the 
cubic  spline  interpolant  should  shed  light  on  the  validity 
of  Marsden's  conjecture  that  the  parabolic  spline  offers 
better  approximation  than  the  cubic  spline  when  functions 
of  the  classes  C^-"-^,  C^^^,  and  C^-^^  are  considered. 

We  first  recapitulate  the  properties  of  the  parabolic 
spline.   Let 

f  6  C[a,b]  ,   f(a)  =  f(b)  , 
I  I  f  I  I  =  sup  C  I  f  (x)  I  :  a  <  X  <_  b  } 
such  that  f  is  extended  periodically  with  period  b  -  a. 

A  function  s  (x)  is  defined  to  be  a  periodic  quadratic 
spline  interpolant  associated  with  f  and  a  partition 

f^i^i^o  if 
(5.1.1) 

a)  s(x)  is  a  quadratic  expression  on  each  (Xj^_-,,  x-)  ; 

b)  s(x)  e  C [a,b]  ; 

c)  s(a)   =  s(b)  ,   s' (a)  =  s' (b)  ; 

d)  s(zj^)  =  f(2^)    i  =  1,  2  .  .  ,  k 
where  z^^^   =  (x^^^+x^) /2. 

The  following  theorem  is  due  to  Marsden  [1974]  and 
was   given   in   Chapter    One    as    Theorem    1.13. 
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Theorem  Let  fx^^l^^g  ^^  ^  partition  of  [a,  b],  f(x)  be  a 
continuous  function  of  period  b  -  a,  and  s  (x)  be  the 
periodic  quadratic  spline  interpolant  associated  with  f 
and  Cx^}^^Q.   Then 

(5.1.2)  lis. II  <  2  ||f||  ,  llsjl  <  2  ||f|i  ; 
I  |e^|  1  <  2   (^(f,h/2)  ; 

I  |e  II  <  3   a)(f ,h/2)  . 
(where  s^   =   s(x^)  and  ej^  =  y(Xj|^)  -  s(Xj^)  ). 
The  constant  2  which  appears  in  the  first  of  the  above 
equations  can  not,  in  general,  be  decreased. 

For  continuous  functions  to  be  "well-approximated"  by 
the  spline  s.  Equations  (5.1.2)  show  that  the  only 
requirement  for  the  partition  is  that  the  length  h  of  the 
largest  subinterval  be  small  enough  that  the  modulus  of 
continuity  of  f  be  small. 

Concerning  s,  we  can  prove  the  following  results. 
These  are  analogous  to  the  results  of  Marsden  given  above 
as  Theorems  1.14  to  1.16  and  improve  upon  the  bounds  he 
derived. 

Theorem  5.1  Let  f  and  f  be  continuous  functions  of  period 
b  -  a.   Then 

(5.1.3)  I  |e(x)  I  I  <  Cq^i  h  |  |f'l  I  ' 

where  ag  =  2/3  -  /13/6  and  Cq  -,  =  1  +  ag  -  Sag  +  4aQ  or 
Cq  -j^  is  approximately  1.0323.  The  analogous  constant  from 
Marsden  was  5/4. 
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Theorem  5.2   Let  f,  f  ,  and  f"  be  continuous  functions  of 
period  b  -  a.   Then 

(5.1.4)  I  |e|  I  <  (1/6)  h2  I  |f  •  I  I  , 

(5.1.5)  I  le^'l  I  <  (9/16)  h  Mf  'I  I  ' 

(5.1.6)  He'll  <  (17/16)  h  I  |f "  I  I  . 

(Marsden's  constant  for  (5.1.4)  was  5/8,  while  in  (5.1.6) 
the  value  was  2)  . 

If  we  make  the  additional  assumption  that  the 
partition  consists  of  equally  spaced  intervals,  then  we 
can  improve  (5.1.6)  to 

(5.1.7)  ||e'  I  I  <  .7431  h  |  |f '  '  |  I  • 

Theorem  5.3  Let  f,  f,  f",  and  f"  be  continuous 
functions  of  period  b  -  a.   Then 

(5.1.9)  I  |e^|  I  <  (1/24)  h^  I  |f ■  "  II  ' 

(5.1.10)  I  le^' I  I  <  (1/6)  h2  I  |f '  ■  •  I  I  , 

(5.1.11)  I  |e|  I  <  (1/24)  h^  Mf "  '  M  ' 

(5.1.12)  He'll  <  (7/24)  h^  |  |f '  '  H  I  > 
(5.1.13) 

I  |e' 

Marsden's  analogous  constants  for  (5.1.9)  to  (5.1.11) 
are  1/8,  1/3,  17/96,  and  11/24  respectively. 

Furthermore,  (5.1.9)  and  (5.1.11)  are  best  possible. 
In  fact  we  also  have  the  exact  pointwise  bound 
(5.1.14)   |e(x)  I  <  |E3(t)  I  h^  I  |f  •  "  II,  x^  <  X  <  x^^;^  , 
where  t  =  (x  -  Xj^)/(x^^-i  -  x^^)  and 


<  [h^/2  +  (h3/3hi2)]  I  |f,,,|  I  ^  ^^__^. 


Q3  (t)  =  1/24  -  t2/4  +  t^/6 
is  the  "Euler  spline"  of  degree  3. 
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The  technique  used  here  is  the  same  as  that  used  in 
the  last  chapter.   For  a  given  partition  subinterval 

f^i'  ^i+l] '  ^^  write 

(5.1.15)   If^^Mx)  -  s^^Mx)|  <  If^^Mx)  -  L^^Nx)! 

+  iL^i) (X)  -  s(i) (X) I 
•where  L  is  a  polynomial  interpolation  of  f.   We  then 
proceed  by  obtaining  pointwise  estimates  of  the  quantities 
on  the  right  hand  side  of  (5.1.15). 

Proof  of  Theorem  5.1 

Given  that  f  and  f  are  continuously  dif f erentiable 
of  period  b  -  a,  we  will  establish  the  following  pointwise 
bound  for  the  parabolic  continuously  dif f erentiable  spline 
s  interpolating  function  values  at  subinterval  midpoints 
z^.      Let  Xj^  _<  X  <_  Xj^^^  and  let 

t  =  (x  -  Xj^)/(x^^^  -  x^)    and  h  =  max{hj_]. 
Then  for  x^^  _<  x  _<  ^i  +  i,    we  have 

(5.2.1)  |f(x)    -   s(x)|    <h    [1    +    t   -    8t2    +    4t2]       llf'll     • 
For    Zj_  +  ]^    1^1    ^i  +  1    i^eplace     t     in     (5.2.1)     by    1     -    t. 
Equation    (5.1.3)    follows   from    (5.2.1). 

In   order    to    establish    (5.2.1)    we    write    for 
Xj^    <^  X   <^  ^i  +  1 ' 

f(x)    -    s(x)     =    f(x)    -   L(x)    +    L(x)    -    s(x) 
where   L(x)    is    the    parabola    matching    f{x^),    f(z.^-|^),    and 
f  (Xj_^-[.)  •      Then 

(5.2.2)  |f(x)    -    s(x)|    <    lf(x)    -   L(x)|    +    |L(x)    -    s (x) I     . 
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We  can  represent  L(x)  as 

(5.2.3)  L(x)  =  f(x^)AQ(t)  +  f  (z^^;L)Aj^(t)  +  f(x^^^)A2(t) 
where  t  =  (x  -  x.)/h-  and 

Ao(t)  =  2  (1/2  -  t) (1  -  t)  , 

A;L(t)  =  4t  (1  -  t)  , 

A2(t)  =  2t  (t  -  1/2)  . 
As  L  reproduces  parabolas  exactly  and  as  the  restriction 
of  s(x)  to  [x^,   x^_^-^]    is  a  parabola,  for  x^^  <_  x  £  ^i  +  i  we 
have 

(5.2.4)  s(x)  =  s(x^)AQ(t)  +  s(z^^^)A]L(t)  +  s  (x^^^^^  ^2  (^)  ' 
As  f(Zj_^-]^)  =  s(Zj^^-|_),  we  have 

(5.2.5)  |L(x)  -  s(x)|  <  |f(x^)  -  s(x^)|  |AQ(t)| 

+  If(x.^l)  -  s(x.^i) I  |A2(t) I 

<  [  |Ao(t)  I  +  |A2(t)  I  }  llejl 

<  |1  -  2t|  Ijeil) 
where  ||ej_||  =  max^<  j_^j^  |  f  (Xj, )  -  s(x^)|  . 

We  have  shown  that 

(5.2.6)  |f (X)  -  s(x) I  < 

|f (X)  -  L(x) I  +  |1  -  2t|  I |e^| I  . 
It  remains  to  bound  |f (x)  -  L(x)  |  and  |  \e^\   |  in  terms  of 

llf'll. 

Marsden  showed  that 

(5.2.7)  I le^l I  <  h  llf'll 

where  h  is  the  maximum  length  of  a  subinterval. 

In  order  to  bound  |f(x)  -  L(x)|  we  resort  to  the 
Peano  theorem.   Defining  g(t)  =:  f(Xj^+hj^t)  =  f(x),  we  have 
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1 

(5.2.8)  f(x)  -  L(x)  =  /  K.  (t,z)g'  (z)dz 

0 
where 

K3^(t,z)  =   (t-2)$  -  Ao(t)  [0-z]0 

-  A-L(t)  [1/2  -  2]° 

-  A2(t)  [1  -  zlj 

and 

1   for  t  >  z 
(t-z)O  = 

0  for  t  <  z  . 

In  order  to  verify  (5.2.8)  ,  one  need  only  expand  the  right 

hand  side  and  integrate  by  parts.   For  0  <.  t  <_  1/2, 

Ki  (t,z)  may  be  written  in  the  more  convenient  form 

K^(t,z)  =  -A-L(t)  -  A2(t)  +  1     for  0  <  z  <  t 

=  -Aj_(t)  -  A2  (t)         for  t  £  z  <  1/2 

=  -A2(t)  for  1/2  1  z  <  1  , 

From  Equation  (5.2.8)  it  follows  that 

1 

(5.2.9)  |f(x)  -L(x)|  <_     j      |K.  (t,z)  |dz  maXQ<^<^|g'(t)  I 

0  

1 

<  h^  J  |K^(t,z)  |dz  maxj^  <j^<j^    |f'(x)| 

0  i i  +  1 

1 

<  h.  /  |Ki(t,z)  |dz  Kf'li  . 

0 

Evaluating  the  integral  in  (5.2.9),  we  have 

1 

(5.2.10)  /  |K, (t,z)|dz  =    t   [1  -  A. (t)  -  A2(t)] 

0 

+  (1/2  -  t)  [A]_(t)  +  A2(t)] 

+  (1  -  1/2)  [-A2(t)] 

=   3t  -  8t2  +  4t^  . 
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Combining  Equations  (5.2.7)  -  (5.2.10)  we  have  for 
0  <  t  <  1/2, 
(5.2.11)    |f(x)  -  s(x)|   <_  [h  (1  -  2t)  + 

hj_(3t  -  8t2  +  4t^)  3   I  |f'  I  I 
<  h  [1  +  t  -  8t2  +  4t^]  I  |f'  I  1 
which  is  precisely  the  desired  result.   The  maximum  of  the 
right  hand  side  of  (5.2.11)  occurs  for  aQ  =  2/3  -   13/6. 
Evaluating  gives  the  value  Cq  -j^. 

Proof  of  Theorem  5.2 
Let  f  be  twice  continuously  dif ferentiable  of  period 
b  -  a  and  let  a  partition 

a  =  Xq  <  z-]_  <  x-j^  <  .  .  Xj^  <  Z-L  +  i  <  Xj^^j^  <.  .  <  Xj^  =  b 
be  given  (where  z^^-^   =    {x^   +  Xj_^-|^)/2,  every  i).   Let  s  be 
continuously  dif f erent iab 1 e  and  a  parabola  on  each 
interval  [x^^,  ^i  +  i^   such  that 

s(Zj^^j^)  =  f(z^^j_),  s(a)  =  s(b),  and  s'(a)  =  s'(b). 
Letting  t  =  (x  -  x^) /h^,    we  show  that 
(5.3.1)    |f (X)  -  s(x)  I  <  Co^2(^)  I  |f "  I  I 

Cq  2(t)  =  h^  C(l  -  2t)/6  +  [  t/(3  -  2t)  -  t^]] 


and  for  z^^^^.  i.  ^  i.  ^i  +  l 

'^0,2(t)     =    Co^2(l-t)     • 
Furthermore    the    maximum    of    Cq   2(t)    is    1/6    and   occurs 

for   t   =    0    and    1. 

As  in  the  proof  of  Theorem  5.1  we  fix  i  and  let  L(x) 

be  the  parabola  satisfying 
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L(Xi)=f(x^)  ,     L(z^^^)=  f(Zi+i)  / 

Then,  proceeding  in  the  same  way  as  before, 
(5.3.2)    |f(x)  -  s(x)|  <  |f(x)  -  L(x)|  +  1 


eill     |1   -    2tl 


We    must    bound     |f(x)    -    L(x)|     and     ||e^||.       We    first 
bound    ||e^||.      From    Marsden    [1974],    we   have 

(5.3.3)  h^    Si_i    +    3(hi    +   hi_i)    s^    +   hi_i    s^+i 

=  4  [h^  f(Zj^)  +  h^_-^   f(Zi  +  i)]  • 
Denoting  f^  =  f(x^)  and  e^^  =  f ^^  -  s^,  we  obtain  from 
Equation  (5.3.3) 

(5.3.4)  h.  e^.i  +  3(hi  +  h^.i)  e^  +  h^.^  e^^^ 

=  h^  fi_i  -  4  h.  f(Zi)  +  3(hi  +  hi_i)  fi 

-  4  hi_i  f(Zi+i)  +  h^-l  fi+1 
=  :  B(f)  . 
As  B  is  identically  zero  for  any  linear  function  f, 
we  have  by  the  Peano  Theorem: 

(5.3.5)  B(f)  =    /  ^^■^K(y)  f"(y)  dy  /  1! 

^i-1 
where 

K(y)  =  Bj^[(x-y)^] 

=  h._i  (Xi^3_  -  y)+  -  4hi.i  (z^^i  -  y)  + 

+  3(hj^  +  h^.^)  (x^  -  y)+  -  4h^  (Zi-Y)  + 
+  hi  (Xi_i  -  y)+ 

and  X  -  y        for  x  >  y 

(X  -  y)+  = 

0  for  X  <  y  . 
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In    order    to    illustrate    the    symmetry   of    the    kernel 
K(y)    about   Xj_,    we   expand    in   terms   of   y   -   x^   to   obtain 

K(y)    =  h^.i    f^i   ~    (y-^i)] 

for   hj^/2   <   y-Xj_   <_  h^^ 

=   h^_;L    t3(y-Xi)    -   h^] 

for    0    <^  y-Xj^    <   hj^/2 
=   h^    [-3(y-Xi)    -   h^.^] 

for   -h^_i/2   <   y-xj_   <    0 
=  h^    [h^_i   +    (y-x^)] 

for   -h^_-]_  1  Y   -   ^i  1  ~hi-i/2 
where  h^^   =  Xj_^-]_  -    Xj_  and  h^_-^  "^   ^i   ~    ^i-l*         ^^   ^^   easily 
seen,    the   sign  of  K(y)    changes  at  y  =  Xj_   +  hj_/3  and 
Xi     -    hi_3_/3. 

From     (5.3.5),     it    follows    that 
(5.3.6)  |h^    e^_i    +    3(hi    +    hi_i)    e.     +    h^.^^    ^i  +  li 


X 


i  +  1 


<  /  |K(y)  1    dy      I  |f"|  I 
^i-1 

<  (h.    +  hi_i)    hi  h^_^    I  |fM/3    . 

I    =    I  |ej_|  I .       Then 

(5.3.7)  I  le^l  I  <  (1/6)  h2  |  |f"  |  I  / 
which  is  the  desired  bound  on  |  je^] ]. 

We  next  bound  |f  (x)  -  L(x)  |  where  L  is  the  parabola 
matching  f  at  Xj_,  z^  +  i/  and  ^i  +  i'  L  can  be  uniquely 
expressed  as 

(5.3.8)  L(x)  =  f(x^)    AQ(t)  +  f(Zi+i)  A^(t) 

+  f(Xi^^)  A2(t) 
where 
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Ao(t)     =    2     (1/2    -    t) (1    -    t)     , 
A^(t)     =    4t    (1    -    t)     , 
A2(t)     =    2t    (t    -    1/2)     . 
Then,  defining  g(t)    =:   f  (x^^  +  hj^t)    =   f  (x)    we  have 

1 

(5.3.9)  f(x)    -    L(x)     =    /    K2(t,z)    g''(z)dz    ,    t    =     (x-x^) /h^ 

0 

where 

K2(t,z)    =    (t   -    z)^    -   AQ(t)[0-z]^ 

-  A^(t)  [1/2  -  z]^  -  A2(t)  [1  -  z]^  . 
Equation  (5.3.9)  can  be  verified  by  integrating  by  parts 
to  obtain  (5.2.8).   For  0  <_  t  <_  1/2,  K2  takes  the  form 

(5.3.10)  K2(t,z)  =  z  (2t  -  1)(1  -  t)  for  t  >  z,  t  <  1/2, 

-  t  [1  +  z(2t  -  3)]     for  t  <  z  £  1/2 

-  t  (2t  -  1) (1  -  z)  for  z  >  1/2,  t  <  1/2  . 
From  (5.3.9),  it  follows  that  for  0  <  t  <  1/2 


(5.3.11)   |f(x)  -  L(x)  |/|  |g"  (t)  I  I  <   /  |K2(t,z)|dz 

0 


-   /  z(2t  -1) (1  -  t)  dz 
0 

l/(3-2t) 
-   /         -t[l  +  z(2t  -  3)]  dz 

X 

1/2 
+   /         -t[l  +  z(2t  -  3) ]  dz 
1/ (3-2t) 


/   2t(t  -  1/2) (1  -  z)dz 
1/2 

-t2  +  [t/(3-2t)]  . 
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Therefore  ,  if  0  <  t  <_  1/2,  we  have 

(5.3.12)  if(x)  -  L(x)  i  <  hi2  [-t2  +  [ t/ (3-2t) ] } 1 | f ' ' (x) || 
We  can  now  assemble  the  parts  to  get  the  pointwise  bound 
(5.3.1).  Using  the  bound  for  1  f  (x)  -  L(x)|  of  (5.3.11) 
and  the  bound  of  (5.3.7)  for  |  |  e^^  |  |  in  the  formula 

|f(x)  -  s(x)|  <  |f(x)  -  L(x)|  +  lleill  |1  -  2t|  , 
we  then  have  for  0  <_  t  _<  1/2, 

(5.3.13)  |f(x)  -  s(x)|  <  [  hi^  [.^2  +  t/(3-2t)] 

+  (1  -  2t)  h2/6  3  I  |f"  1  1 
which,  as  h^  <  h,  immediately  implies  (5.3.1).  The  result 
for  1/2  1  t  <_  1  follows  by  symmetry.  It  remains  only  to 
be  shown  that  the  maximum  of 

is  h^/e  and  occurs  for  t  =  0.  To  see  this,  expand  CQ^2(t) 
at  0  as 

^0,2^^^  =  ^0,2^°)  ^  ^  ^0,2'(0^  ^  <^^/2)  Co^2"<y) 
where  0  <  t  <  1/2  and   0  <  y  <  t.   It  is  not  hard  to 
verify  that  the  last  two  terms  of  the  above  expression  are 
negative,  and  hence  the  maximum  occurs  at  t  =  0.   This 
completes  the  proof  of  Equation  (5.1.4). 

We  next  show  Equation  5.1.5.   From  Marsden,  we  have 
the  tridiagonal  system  matching  spline  derivatives, 
(5.3.14)   h._i  s._i'  +  3(h.  +  h._^)  s^'  +  h.  s^^i' 

=  8  [f (z^+^)  -  f (z^)] 
or  equivalently , 
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(5.3.15)  h^_i   e^_^'    +    3(h^    +   h^.^L^    ^i '    ^   ^i   ^i  +  l' 

=   hj__^    f^.l'    +    8f(Zi)    +    3(h^    +   h^_i)    f^' 

-    8f(Zi^i)    +   hi    fi^i' 
=:    B^lf)     . 

where 

Bi[(x-y) J    =   h^  z^+l   <   y   <   x^^i 

=    8(y-X3^)  -    3hi  ^i    <    Y   <    ^i  +  i 

=    8(y-Xj^)  +    3hj^_3_        Zi    <   y   <   Xi 

=   -hi_i  H-1   <y  1   ^i 

and  the  last  equality  of  (5.3.15)  can  be  verified  by 
integrating  by  parts  and  using  the  fact  that  B -^  is 
identically  zero  if  f  is  a  linear  function.  From  (5.3.15) 
follows 

(5.3.16)  |hi_i    e^.^'    +    3(hi    +    h^_^)    e^'    +    h^    e^^^'l 

<  /''^'■^lBi[(x-y)  J  |dy    1  |fl  I 
^i-1 

<  (9/8)     (h^^    ^    h^_^^)     \\f"\\    , 

and   hence 

(5.3.17)  I  le^'  I  I    <     (9/16)    h    Mf '  '  |  I     ' 
which   is   precisely    (5.1.5)  . 

We    next   establish   a   pointwise   bound   on    |e'(t)|.      We 
show  that  for  Xj^   <_  x  <_  Xj_^-|_,   and  t   =    (x  -    Xj_)/hj_, 

(5.3.18)  |f'(x)    -    s'(x)|     <    C(9/16)h 

+    2    h^    t(l    -    t)  }     I  |f"|  1     • 
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If  we  maximize  the  right  hand  side  of  (5.3.18)  by  taking 
h.  =  h  and  t  =  1/2,  then  (5.1.6)  is  an  immediate 
consequence  of  (5.3.18). 

To  establish  (5.3.18),  we  let  J(x)  be  the  unique 
parabola  satisfying 

(5.3.19)  J(z^  +  ]_)   =  f(z^  +  ^)  ,   J'(Xj_)  =  f'(x^)  , 
J'(Xi^l)  =   f  (Xi^i)  . 

Then 

J'  (X)  -  s-  (X)  =  (1  -  t)  ej_'  +   t  e^^^L* 
and 

(5.3.20)  |f '  (X)  -  s'  (X)  I 

<  |f '  (X)  -  J'  (X)  I  +  |J'  (X)  -  s'  (X)  I 

<  |f '  (X)  -  J'  (X)  I  +  I  le^'  I  I  (|l-t|  +  |t|3 

<  |f'  (X)  -  J'  (X)  I  +  I  le^'  I  I  . 

As  we  already  have  an  estimate  for  |  |  e^^'  |  |  ,  we  need 
only  estimate  |f'(x)  -  J' (x)  |  .  It  is  easy  to  see  by 
integration  by  parts  that 

1 

(5.3.21)  f'(x)  -  J' (x)  =   /  K(t,z)g' ' (z)dz  /  h^ 

0 

where  g(t)  =:  f(Xj_+hj^t)  =  f  (x)  and  where 
K(t,z)   =     1  -  t    for  t  >  z 

-  t     for  t  <  z  , 
From  (5.3.21)  follows 

1 

(5.3.22)  If'(x)  -  s'(x)|  <   /  |K(z,t)|dz  ||g"ll  /  ^i 

0 

and  evaluation  of  (5.3.22)  gives 
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(5.3.23)  |f'  (X)    -    J'(x)|    <      h^    2t    (1    -    t)     ||f"|l     ■ 
Applying    (5.3.23)    and    (5.3.17)    in    (5.3.20)    gives 

(5.3.24)  |f ' (X)    -   s' (X) I    < 

Ch^    2t    (1    -   t)    +       (9/16)    h    3       I  |f"  I  I     ' 
which   is    the   desired   pointwise   bound. 

In  order  to  improve  the  bound  for  the  case  of  evenly 
spaced    knots     we    return    to    the    use    of    the     parabola     L 
satisfying 

L(x^)    =   f(x^),    L(z^^-^)     =    f(Zi  +  i)     , 

L(Xi  +  i)    =    f(Xi  +  i)     . 

Then  for  arbitrary  i  and  Xj_  <_  x  _<  x^^-j.'  ^^  have 

(5.3.25)  |f ' (X)  -  s' (x) I 

<  |f ' (X)  -  L' (X) I  +  |L' (X)  -  s' (X) I 

<  |f •  (X)  -  L'  (X)  I 

+  (1/h^)  le^  Ao'(t)  +  e^^i  A2'(t)| 

<  |f •  (X)  -  L-  (X)  I 

+  (| le^l |/h^)  {|Ao' (t) I  +  IA2' (t) |} 
where  t  =  (x-Xj_)/hj_.  Differentiating  equations  (5.3.8) 
gives 

(5.3.26)  IAq' (t) I  -  |3  -  4t|  , 
IA2' (t)  i  =  |1  -  4t|  . 

Recalling    from    Equation    (5.3.7)     that 

lie.  II    <     (1/6)    h2    llf'M     . 
we   have 

(5.3.27)  |f(x)    -    s'(x)|    <    |f'(x)    -   L'(x)| 

+     (h2/6h^)     { |Aq'  (t)  I  +  IA2'  (t)  |3     ||f  "  I  I 
where 
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|Aq' (t) I  +  IA2' (t) I  =  4  -  8t  ,   0  <  t  <  1/4 

=  2       ,   1/4  <_  t  £  3/4 
=  8t  -  4  ,   3/4  <  t  <  1  . 
To  obtain  a  pointwise  bound  we  need  only  bound 
|f' (X)  -  L'(x)|.   Differentiating  (5.3.9),  we  obtain 
(5.3.28)   If  (X)  -  L'  (X)  I  < 


h.   /  |K  <1'°)  (t,z)  |dz  I  \f"   (X) 
^   0 


where  for  0  <  t  <  1/2 


K2^^'°^  (t,z)  =  z(3  -  4t)  ,       t  >  z 

-  1  +  z(3  -  4t)  ,  t  <  z  <  1/2 

(1  -  4t) (1  -  z)  ,  1/2  <  z  <  1  . 

The  only  sign  changes  occur  for  z  =  1/(3  -  4t)  and 

0  <_  t  <  1/4   and  along  the  lines  t  =  1/4  ,  z  <_  1/2,  and 

z  =  t.   Evaluating  the  the  integral  of  (5.3.28)  for 

0  <_  t  <_  1/4  then  gives 

1  t 

(5.3.29)    /  |K2^-'-'°Mt,z)  |dt  =    /  z(3  -  4t)dz 
0  0 

l/(3-4t) 
-  /         [-1+ z(3-4t)]dz 
t 

1/2 
+  /        [-l  +  z(3-4t)]dz 
l/(3-4t) 

1 
+  J    (1  -  4t)  (1  -  z)dz 
1/2 

=  t2(3  -  4t)  -  2t  +  1/ (3  -  4t)  . 

For  the  interval  1/4  <_  t  <_  1/2,  the  only  sign  change 

is  the  line  2  =  t.  We  obtain 
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1  X 

(5.3.30)  /  |K,(1'°) (t,z) |dt  =    /  z(3  -  4t)dz 
0  0 

1/2 

-  /    [-1  +  z(3-4t) ]dz 

X 

1 

-  /  (1  -  4t) (1  -  z)dz 

1/2 

=  t^O  -  4t)  . 
The  expression  of  (5.3.30)  is  monotone  increasing  for 
1/4  1  t  _<  1/2.   We  have  shown  that 

(5.3.31)  If  (X)  -  L'  (X)  |/|  |f  •  I  I 

h^    Ct2(3  -  4t)  -  2t  +  l/(3-4t)}  , 

0  <  t  <  1/4 

<  2 

h^  ^^(3  -  4t)  , 

1/4  <  t  <  1/2 
Combining  (5.3.27)  and  (5.3.31)  we  have 

(5.3.32)  |f'  (X)  -  s'  (X)  |/|  |f  "  I  I  < 

hj^[t2(3-4t)  -  2t  +  l/(3-4t)] 


+  (h2/3hi) (2-4t)  , 


0  <  t  <  1/4 


<  hj_[t2(3-4t)]  +  h2/3h^ 


1/4  <  t  £  1/2  . 
As  usual  these  results  can  be  extended  by  symmetry  to 
the  interval  1/2  <.  t  <_  1.  The  bound  given  in  (5.3.32)  is 
monotone  decreasing  from  0  to  1/4  and  increasing  from  1/4 
to  1/2.  For  an  equally  spaced  partition,  and  t  near  1/2, 
(5.3.32)  is  quite  a  bit  smaller  than  (5.3.23).  For 
instance  for  t  =  1/2,  we  have  1/3  +  1/4  versus  17/16  from 
(5.3.23).   If  we  rewrite  (5.3.23)  as 


98 


(5.3.33)  |f'(x)  -  s'(x)|  <  C2  i(t)     ||f' 'II 
where 

C2^l(t)  =  h^  2t  (1  -  t)  +  (9/16)  h 
and  write  (5.3.32)  as 

(5.3.34)  |f'(x)  -  s'(x)|  <C2^i(t)  ||f' '11 

where  C2  i(t)  is  the  right  hand  side  of  (5.3.32),  then  we 
have  also  the  pointwise  bound 
(5.3.34)   |f • (X)  -  s' (X) I  < 

min[C2^l(t)  ,  C2^i(t) }  |  |f  "  |l  • 
For  h  =  h- ,  the  maximum  of  the  right  hand  side  of  (5.3.34) 
is  approximately  .7431  h  |  |f''|  1  and  occurs  for  t 
approximately  equal  to  .10038. 

This  completes  the  proof  of  Theorem  5.2. 

Proof  of  Theorem  5.3 
We  first  show  that  if  f,  f,  f",     and  f"  are 
continuous  and  of  period  b  -  a,  then 

(5.4.1)  ||f  -  s|  I  <  (1/24)  h^  I  |f"  ■  II 

where  s  is  the  parabolic  and  periodic  (once  differen- 
tiable)  spline  interpolating  subinterval  midpoints. 
Furthermore,    "1/24"    cannot   be    improved. 

We    proceed    by    demonstrating    (5.1.14).       For    a    given 
partition   and   subinterval    [x^,    x.j_  +  i]  r   we   write 

(5.4.2)  |f(x)    -    s(x)|    <     |f(x)    -   L(x)|    +     |L(x)    -    s(x)| 
where   L(x)    is    the   unique    parabola    satisfying 

(5.4.3)  L(x^)     =    f(x^)     ,    L(z^^;i^)     =    f(Zi  +  i)     r 
L(^i+l)     =    f(^i+l) 
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where     z  ^  +  ^     =     (^i  +  l     ^    Xj_)/2.       L(x)     may    be    uniquely 
expressed   as 

(5.4.4)  L(x)     =    f(x^)    Ao(t)     +    f(Zi+i)    A^(t) 

+    f(xj_^j_)    A2(t) 
where   t   =    (x   -   x^) /h^,    and 

AQ(t)  =  2(1/2  -  t) (1  -  t)  , 
A-L(t)  =  4t(l  -  t)  , 
A2(t)  =  2t(l/2  -  t)  . 
Proceeding  by  using  the  Cauchy  formula  one  obtains 

(5.4.5)  |f (x)  -  L(x) I  < 

(1/6)  |t(l/2  -  t)  (1  -  t)  I  I  |f'  "  I  I  . 
In  order  to  bound 

(5.4.6)  |L(x)  -  s(x)|  =  I  [f(x^)  -  s(x^)]    AQ(t) 

+  [f  (x^^^)  -  s(Xi^.^)]  A2(t)  I 
<  lleill  f  |Ao(t)|  +  |A2(t)|  3 
=  lleill  |1  -  2t|  , 
we  must  bound  ||ej_||  =  max^_-j^  2   k'-l^il-'  ^l^^re 
e^  =  f(x^)  -  s(x^).   To  bound  |  |ej_|  |,  we  resort  in  turn  to 
the  tridiagonal  system  of  (5.3.4)  , 

(5.4.7)  h^  e^.i  +  3(hi  +  h^.^)  e^  +  h^_^  e^+^ 

=  h^f^.i  -  4hif(2i)  +  3(hi  +  hi_i)fi 

-    4^1-1^(^1+1)     ^    h^.^f^^i 
=:    Bi(f)     . 
B2_(f)  r    SO    defined,    is    a    linear   functional    identically 
zero  for  polynomials  of  degree  two  or  less.     We  thus  have 

(5.4.8)  Bi(f)     =      /    ^^^K(y)f "  ■  (y)dy    /    2! 

^i-1 
where 


100 


K(y)     =    B^[(x-y)^2] 

=   h^_i    [x^+i   -  y]2      -   4hi_i    [    z^^i   -  y]2 
+    3(hi    +   h^_^}     [x-    -   y]^2 

-  4hi    [Zi   -   y]2    +   h^.i    [x._i   -   y]2 

=   h^_^    [h^   -    (y-x^)]2    ,      hj_/2   <   y-x^   <   h^ 
^i-l^^i   "    (y-^i']^   -   4h^_^[hi/2   -    (y-Xi)]2    , 

0    <   y-Xj_    <      h^/2 
~^i[^i-l    ■"    (y-^i)^^    +   4hi[hi_^/2    +    (y-Xi)]2    , 

-hj__^/2    <   y-x^    <    0 
~^ithi-i    +    (y-x^)]2    ,    -h^_^<   y-x^    <    -h^_i/2    . 
Then 

(5.4.9)  Ih^e^.i    +    (h^    +   hi_l)ei    +   h^.^e^^^Ll 

<      /    ^^^iK(y)  |dy    I  |f' "  M    /    2! 

-  (1/12)     (hi_ihi^    +   h^h^.i^)     ilf  "I  I     . 
If   we   take    i    so   that    \e\\    is   maximal,    then  we   have 

(5.4.10)  2(h^    +   hj__^)     I  le^l  I    < 

=    (1/12)     (h^.^hi^    +   h^h^.i^)     I  |f  "I  I 
and   hence 

(5.4.11)  ||e.|| 

<  (1/24)     (h^_^h^^    +   h^h^_^^)/(hi+hi_^)     l|f' "II 

<  (1/24)    h^     I  |f "  ■  I  I     ■ 
Combining    (5.4.6)    and    (5.4.11),   we  have 

(5.4.12)  |f(x)    -    s(x)|    <     lf(x)    -   L(x)|    +    11^^11     |1    -2t| 
<     [|t(l/2-t)  (1-t)  |/6     +     (1/24)  |l-2t|3  h^  1  If  "I  I 

=    lQ3(t)  I    h^     I  |f  "  '  I  I 
where    Q3(t)     is    the    Euler    spline    of    degree    three.    Equation 
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(5.4.12)  is  precisely  (5.1.14).  From  it  follow  also 
(5.1.9)  and  (5.1.11). 

To  see  that  (5.1.14)  cannot  be  improved,  consider  the 
"Euler  spline"  Qn(x)  constructed  by  integrating  a  constant 
n  times  so  that  the  nth  integration  is  odd  for  n  odd  and 
even  for  n  even.  On  the  unit  interval,  the  first  few 
Euler  splines  are 

(5.4.13)  Qq{x)    =    1    , 
Q^(x)  =  X  -  1/2  , 
Q2(x)  =  x^/2    -   x/2  , 

Q3 (X)  =  x^/6  -  x2/4  +  1/24  . 
These  can  be  compared  to  the  even  Euler  polynomials  given 
in  Chapter  Two. 

If  we  extend  the  Euler  splines  to  the  real  line  by 
setting 

(5.4.14)  Q^(x)    =    (-1)  JQj^(t  +  j) 

0  <_  t  <_  1  and  j  integer,  then  Q  (x)  is  n  -  1  times 
continuously  di f f eren t iab 1 e  and  piecewise  n  times 
continuously  dif f erentiable.  Q  is  of  period  2  and  has 
nth  derivative  of  plus  or  minus  one.  Q  is  thus  a  member 
of  the  class  of  functions  with  nth  derivative  piecewise 
dif f erentiable.  The  third  derivative  of  Q3  can  be  repre- 
sented as  the  pointwise  limit  of  the  the  third  derivative 
of  a  sequence  ff^}  of  three  times  continuously  differen- 
tiable  functions  which  converge  uniformly  to  Qo.  Further- 
more each  of  the  f ^^  have  third  derivative  bounded  in 
absolute    value    by   one. 
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Restricting  Q3 (x)  to  any  interval  [0,  2k],  consider 
the  once  continuously  dif f erentiable  spline  s  parabolic  in 
each  interval  [i,  i+1] ,  and  satisfying 

(5.4.15)  s(z^+i)  =  03(21+1)  =  0  ' 

s(0)  =  s(2k)  ,  s' (0)  =  s' (2k)  . 
The  spline  s  thus  defined  is  identically  zero.   It  is  not 
hard  to  see  that  the  maximum  error  occurs  at  the  integer 
knots  and  is  1/24.   In  fact  we  have  shown  that  Q3  (x)  is 
actually  a  pointwise  exact  bound. 

We  next  demonstrate  Equation  (5.1.10), 
(5.1.10)   I  |e^'  I  I  <  (1/6)  h2  I  |f "  •  I  I  , 

where  ||ej_'|]  =  maxp,,  j_<  j,  |  f '  (Xj_)  -  s'(xj_)|.   The  proof 
uses  the  same  functional  B^^  that  we  used  in  proving 
(5.1.5).   As  in  that  case  we  have 

(5.4.16)  hi.i  ei_i'  .  3(h.  .  hi_i)  e^'  .  h^  e^^^' 

=  h^_l  f^_-^'    +    8f(Zi)  +  3(h^  +  h^.-L)  f^' 
-  8f(z.^i)  +  h.  f.^i' 

=:       Bi(f)     , 
which    is    identically   zero   for   all   polynomials   of  degree 
two   or    less.      Hence 

Xj 

x^ 
where 


(5.4.17)        |B^(f)|     <      /    ^^    \Bj_[{x-Y)^^]\dY    |lf'"l|/2! 

•i-1 


Bi[(x-y)J]  =  2h^(x^^-^  -  y)^  -  8(z^^^  -  y)  ^ 
+  6(hi  +  h^.;^)  (y-x^)^  +  S{Y-z^)l 
+    2h^_^(y-x^_^)^ 


103 


=  2h^2  _  2hj^(y-x^)  , 

h^/2  <  y-x^   <   h^ 
6h^(y-X3^)  -  8(y-x^)2  , 

0  <  y-x^  <  hj^/2 
-eh^_-^(Y-x^)    -    8(y-x^)2  , 

Conveniently,  the  above  kernel  is  positive.  Evaluation  of 
the  integral  in  (5.4.17)  is  thus  straightforward,  leading 
to 

(5.4.18)  Ih^.^L  ^i-l'  +  3(h^  +  Yi^_^)    e^'  +  h^  ^i  +  i    \ 

<  (1/3)  (h^^  +  h^.i^)  I  |f  •■•  I  I  . 

Assuming  that  i  is  such  that  |ej^'|  attains  its  maximum  we 

then  have 

3       3 
(5.1.10)   I  |e-  '  I  I  <  (^i   -^  ^i-1   )   ||f"'ll 
6  (h^  +  h^_^) 

<  (1/6)  h^  I  |f • • •  I  I  . 

In  order  to  extend  the  bound  (5.1.10)  to  the  entire 
interval,  choose  any  subinterval  [x^,  ^^+i]  of  the  given 
partition  and  consider  the  line  J'  interpolating  fj_'  and 
f^  +  j^'.   J'  may  be  represented  as 

(5.4.19)  J'(x)  =  (1  -  t)  f^'  +  t  f^^3^'  . 
By  the  triangle  inequality,  we  have 

(5.4.20)  |f' (x)  -  s'(x)|  <  |f'(x)  -  J'(x)| 

+    |J'  (X)    -    s-  (X)  I     . 
As    f    is    twice   continuously  dif f erentiable ,    we   have    the 
well-known  inequality 
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(5.4.21)  |f •  (X)  -  J'  (X)  I  <  h^2  ^(1  _  ^)  I  |f ,  ,  ,  I  I  /  2! 
As  both  J'  and  s'  are  lines  on  [x^,   Xj^^-,  ]  '  ^^  have 

(5.4.22)  |J'(x)  -  s'(x)|  <  [f^'    -   s^']     (1  -  t)  + 


<  llei'll  C  |1  -  t|  +  |t|] 

<  llei'll 

<  (1/6)  h^  I  |f' ' •  I  I  . 

Adding  (5.4.21)  and  (5.4.22)  gives  the  desired  formula 

(5.4.23)  jf  (x)  -  s'(x)|  <  [1/6  +  t(l-t)/2]  h^  ||f'''|| 

<  (7/24)  h^  ||f'  •  •  II  . 

Several  further  refinements  of  this  argument  are  possible. 
This  particular  bound  may  be  worthy  of  further  study  in 
the  future. 

We  next  wish  to  bound  f''(x)  -  s"(x).  Choosing  the 
arbitrary  partition  subinterval  [x-,  Xj^i],  we  consider 
the  parabola  L  matching  f  at  x^,  ^i  +  ir  and  ^i  +  i'  By  the 
triangle  inequality  we  have 

(5.4.24)  I  f "  (X)  -  s"  (X)  I 

<  I  f "  (x)  -L"  (x)  I  +  |L"(x)  -s"(x)| 


<  |f ' ' (x)  -  L' 

+  (l/h^2) 

<  |f • ' (X)  -  L' 

+  (l/h^2) 

<  |f ' • (X)  -  L' 


(X)  I 

e^  Aq'  '  (t)  +  e^^^L  ^2*  '  ^^^  I 

(X)  i 

|eil  I  f|Ao"(t)  j  +  |A2"(t)  |3 


+  (l/h^^)  (h3/24) 


f  •  •  •  I  I  [4  +  4] 
, 3  / TV,  2 ' 


<  |f''(x)  -  L''(x)|  +  (h-^/3h^^)  ||f'''||  . 
The  bound  on  f"(x)  -  L''(x)  is  obtained  by  the  Peano 
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theorem  technique.   We  have  for  0  <_  t  <_  1  and 
f  e  C"  [x^,Xj_^^]  ,  that 
(5.4.25)   |f' ' (X)  -  L' ' (X) I  < 

"^   (2,0) 


h^    |K^^'^'  (t,z)  |dz  Mf  "  I  1/21 
0 


where 


K(t,z)  =   (t  -  z)J  -  Ao(t)  [0  -  z]2 

-  A3_(t)  [1/2  -  z]2  -  A2(t)  [1  -  z]l 
-   (t  -  z)2  -  (1/2  -  z)l    4t(l  -  t) 
-  (1  -  z)2  2t(t  -  1/2) 
and  for  0  <  t  ^  1/2, 

K<2/0) (t,z)/21  =  2z2  t  >  z 

-1  +  2z2         t  <  z,  z  <  1/2 
-  2(1  -  z)2      t  <  z,  z  >  1/2  . 
The  first  of  these  two  terms  is  positive  and  the  last  two 
are  negative.   Evaluation  of  the  integral  of  (5.4.25)  is 
straightforward,  giving  for  0  <_  t  _<  1/2, 

(5.4.26)  If  '  (x)  -  L'  •  (X)  I  < 

h^  i  (4t^/3)  -  t  +  1/2}  I  |f "  '  i  I  • 
Using  (5.4.26)  in  (5.4.24)  gives  for  0  <  t  <_  1/2, 

(5.4.27)  |f''(x)  -  s''(x)|  <  [hj_[(4t^/3)  -  t  +  1/2] 

+  h^/3h^23  I  |f .  .  .  I  I 

<  [hi/2  +  h^/Sh^^^  I  |f •  ■  I  I  , 
which  holds  also  for  1/2  £  t  <_  1. 

Using  the  linear  interpolation  of  f ^'  and  f  j_  +  ]_'  and 
the  usual  triangular  inequality,  we  may  obtain  the 
alternate  estimate 
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(5.4.28).      |f''(x)    -    s''(x)|    <    [hj_[l/2    +    t(l    -    t)  ] 

+    h^/Shj^]     I  If  •  '  '  I  I 
which  when  h  is   larger  than  h^  sometimes  offers   a   lower 
estimate   of   the   error.      As    the   proof    is   very   similar   to 
those   already   given,    we   omit   the   details. 


CHAPTER    SIX 
CONCLUDING    REMARKS 

In  the  above  chapters,  we  have  given  bounds  for  the 
error  of  approximation  of  several  polynomial  and  spline 
interpolations.  Rather  than  restate  the  theorems  proved 
in  previous  chapters,  we  will  try  to  indicate  what  further 
work    is    possible   and   desirable. 

The     work     done     here     has     provided     some     of     the 

motivation     for     the     work     of     Bojanov     and     Varma     (in 

preparation)    extending  Cauchy's   formula  for  the  error  of 

polynomial    interpolation    to    an    expression    for    derivative 

error.       They   proved   the    following   theorem. 

Theorem    6.1    (Bojanov   and   Varma)    Let    f   6    c^"^"^-"-^  [a,b]    and 

let  L[f,x]    be   the   polynomial   of  degree  n  interpolating  f 

at  n  +   1   points.      Then   for   i   <   n   +    1,    we  have 

(6.1.1)  |f <^) (X)    -   L<^^  [f,x] I    < 

II    n.l^'Ml     ||f^^''M|/(n-l)! 

where 

n 
j,,l(x)     =  (X    -    X.  )      . 

1  =  0 

Furthermore,     (6.1.1)     continues    to    hold    for    the    case    of 

Hermite    interpolation     (when       ^^(x)     is    appropriately 

redefined).       In    particular,     let    f    e    c(2"'Mo,l]     and    let 

Vo    _i  by   the   two-point  Hermite   interpolation   satisfying 
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V2^_l(0)  =  f(0)  ,  V2^_i(l)  =  f(l)  , 

for  i=l,  2,  ..,m-l.   Then 
(6.1.2)    |f (^^  (X)  -  V2m-i^^^  (X) | 

<  I  |d_^  x"'(l-x)"'|  I  I  |f  (2m)  I  |/(2ni)!  . 
dx^ 
Equation  (6.1.2)  generalizes  the  results  of  Birkhoff  and 

Priver  [1967].  A  further  generalization  of  (6.1.1)  could 
hold  for  well-posed  two-point  Birkhoff  interpolation  and 
thus  give  the  norm  bounds  of  Chapters  Two  and  Three. 

An  alternate  approach,  suggested  in  conversation  with 
Garrett  Birkhoff,  is  to  compute  pointwise  bounds  for  the 
derivatives  of  polynomial  interpolation  by  the  automated 
approximate  evaluation  of  Peano  kernels.  A  computer 
routine  that  accomplishes  the  necessary  evaluations  is 
given  by  Howell  and  Diaa  (available  on  request).  Similar 
routines  could  automate  some  other  sorts  of  error  bounds. 
Some  of  the  possibilities  are  detailed  in  Howell  and  Diaa. 

We  next  discuss  further  possibilities  for  the  study 
of  spline  error  bounds.  Hopefully,  knowledge  of  the  errors 
of  such  error  bounds  should  aid  in  comparing  various 
spline  operators.  For  instance,  for  the  class  of 
functions  three  times  continuously  dif f erentiable ,  the 
periodic  parabolic  splines  were  shown  in  Theorem  5.3  to 
have  bounds  very  close  to  the  cubic  splines  of  the  type 
discussed  by  Varma  and  Katsifarakis  (Theorem  1.17).  On 
the  other  hand,  many  of  the  parabolic  spline  bounds  can  be 
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bettered  by  a  C^^'  quartic  spline  which  will  be  discussed 
in    later    work. 

Many  other  spline  error  bounds  may  be  amenable  to  the 
same  techniques  employed  here.  Among  these  are  the 
derivative  bounds  for  the  C^  '  quartic  spline  of  Chapter 
Four,  bounds  for  the  lacunary  quintics  discussed  by  Meir 
and  Sharma  [1973],  and  bounds  for  the  local  scheme 
proposed   by   Prasad   and   Varma    [1979]. 

Even  the  extensively  studied  cubic  splines  require 
further  study  along  these  lines.  The  techniques  used  here 
might  yield  good  results  when  the  interpolated  functions 
are  once  or  twice  continuously  dif f erentiable.  Even  the 
second  derivative  bound  given  by  Hall  and  Meyer  [1976]  for 
c'  '  cubic  splines  interpolating  C^  '  functions 
dif f erentiable  is  not  exact.  Finally  of  course  it  would 
be  of  interest  to  develop  a  method  (even  if  merely 
numeric)  of  deriving  optimal  bounds  for  spline 
interpolation  in  each  case  that  the  problem  makes  sense. 
Another  type  of  error  problem  which  is  accessible  by 
techniques  similar  to  the  ones  given  here  is  the  error  of 
the  quadrature  associated  with  any  given  piecewise 
polynomial    interpolation. 
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